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ABSTRACT 

From a sample of 72 galaxies with reliable supermassive black hole masses Mbh, we derive the Mbh~ 
(host spheroid luminosity, L) relation for i) the subsample of 24 core-Sersic galaxies with partially- 
depleted cores, and ii) the remaining subsample of 48 Sersic galaxies. Using if s -band Two Micron 
All Sky Survey data, we find the near-linear relation Mbh oc £^ 10±0 - 20 for the core-Sersic spheroids 
thought to be built in additive dry merger events, while Mbh oc i^ 73±0 - 55 fo r the Sersic spheroids built 
from gas-rich processes. After converting literature B-band disk galaxy magnitudes into inclination- 
and dust-corrected bulge magnitudes, via a useful new equation presented herein, we obtain a similar 
result. Unlike with the Mbh-(velocity dispersion) diagram, which is also updated here using the same 
galaxy sample, it remains unknown whether barred and non-barred Sersic galaxies are offset from 
each other in the M^-L diagram. 

While black hole feedback has typically been invoked to explain what was previously thought to 
be a nearly constant Mbh/Ms p hcroid mass ratio of ~0.2%, we advocate that the near-linear Mbh - 
L and Mbh _ Ms p horoid relations observed at high masses may have instead largely arisen from the 
additive dry merging of galaxies. We argue that feedback results in a dramatically different scaling 
relation, such that black hole mass scales roughly quadratically with the spheroid mass in Sersic 
galaxies. We therefore introduce a revised cold-gas 'quasar' mode feeding equation for semi-analytical 
models to reflect what we dub the quadratic growth of black holes in Sersic galaxies built amidst 
gas-rich processes. Finally, we use our new Sersic M^-L equations to predict the masses of candidate 
intermediate mass black holes in almost 50 low luminosity spheroids containing active galactic nuclei, 
finding many masses between that of stellar mass black holes and supermassive black holes. 
Subject headings: black hole physics — galaxies: evolution — galaxies: nuclei 



1. INTRODUCTION 

The (black hole mass)-(host spheroid dynamical mass) 
Mbh-Ms p h,dyn relation (e.g. Marconi & Hunt 2003; 
Haring & Rix 2004) has recently been revised such 
that core-Sersic galaxies (whose spheroidal components 
are thought to be built in simple additive dry merger 
events) define a near-linear relation while Sersic galax- 
ies (whose spheroidal components are likely built amidst 
gas-rich processes) define a near-quadratic relation (Gra- 
ham 2012a). Graham (2012a, b) additionally predicted 
that the previously established log-linear, i.e. single 
power-law, relation between black hole mass and host 
spheroid luminosity, L, must also require a significant 
modification such that Mbh oc L 10 for the spheroidal 
component of luminous (Mb ^ —20.5 mag) core-Sersic 
galaxies while Afbh oc L 25 for the fainter Sersic spheroids 
— if their dynamical mass-to- light ratio scales with L 1 / 4 . 
If (M/I/)d yn amicai oc L 1 ! 3 then one expects Mbh oc L 8 ^ 3 . 
Here, for the first time, we show that the M^-L relation 
is indeed bent, and we provide new expressions to pre- 
dict black hole masses in core-Sersic galaxies and Sersic 
galaxies using the luminosity of their spheroidal compo- 
nent. 2 

From a literature sample of eight galaxies, Kormendy 

1 Corresponding Author: AGraham@swin.edu.au 

2 Readers unfamiliar with the core-Sersic model (Graham et al. 
2003; Trujillo et al. 2004) or the Sersic (1963) model, for the clas- 
sification of galaxies, may like to refer to the review article by 
Graham & Driver (2005). 



& Richstone (1995, their figure 14) revealed an appar- 
ently linear correlation between black hole mass and the 
brightness / mass of the host spheroid. As the sample 
size grew, researchers began to apply a single log-linear, 
i.e. power-law, relation to the data (e.g. Richstone et 
al. 1998; Kormendy & Gebhardt 2001; McClure & Dun- 
lop 2002). However the samples remained dominated by 
fairly luminous spheroids and the applicability of the new 
M hh -L relations (e.g. McClure & Dunlop 2004; Marconi 
& Hunt 2003; Graham 2007a) to low-luminosity systems 
was not securely established. Indeed, Graham (2007a) 
noted that the M^-L and M^-cr relations (Ferrarese 
& Merritt 2000; Gebhardt et al. 2000) cannot both be 
described by a power-law because the relationship be- 
tween luminosity and velocity dispersion, er, for elliptical 
galaxies is not described by a single power-law. 

As reviewed in Graham (2012c), the elliptical galaxy 
L-a relation has long been known to be 'broken', such 
that L oc er 5 at the luminous-end (Mb < —20.5 mag: 
Schechter 1980; Malumuth & Kirshner 1981; Von Dei- 
Linden et al. 2007; Liu et al. 2008) and possibly even 
as steep as L tx a 6 (Lauer et al. 2007), while L tx a 2 
at intermediate and faint luminosities (e.g. Davies et al. 
1983; Held et al. 1992; de Rijcke et al. 2005; Matkovic 
& Guzman 2005; Kourkchi et al. 2012). Galaxy sam- 
ples containing both bright and intermediate luminosity 
elliptical galaxies will therefore naturally yield L-a re- 
lationships with exponents around 3 to 4 (e.g., Faber 
& Jackson 1976; Tonry 1981; de Vaucouleurs & Olson 
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1982; Desroches et al. 2007). It should be noted that this 
change in slope of the L-a relation at Mg w —20.5 mag 
(Sersic n w 3-4) is not related to (i) the onset of the coex- 
istence of classical bulges and pseudobulges at nbuigc J$ 2, 
nor (ii) the alleged divide at Mb = — 18 mag (n ga i axy < 
1-2) between dwarf and ordinary elliptical galaxies (Kor- 
mendy 1985; Kormendy et al. 2009). Instead, the change 
in slope at Mb ~ —20.5 mag coincides with the division 
between core- Sersic galaxies (which have a central deficit 
of light relative to the inward extrapolation of their outer 
Sersic light profile) and Sersic galaxies 3 which do not 
(Graham & Guzman 2003; Graham et al. 2003; Trujillo 
et al. 2004; Matkovic & Guzman 2005). If (non-barred) 
Sersic and core-Sersic galaxies follow the same Mbh oc a° 
relation (see section 13. ip , they therefore obviously can- 
not follow the same M^-L relation because they do 
not follow the same L-a relation. Simple consistency 
arguments dictate that the core-Sersic galaxies should 
follow the relation Mbh oc L 10 while the (non-barred) 
Sersic galaxies should follow the relation Mbh oc L 25 (or 
M bh oc L 3 if M bh oc a 6 ). 

While the core-Sersic model has provided a new means 
to identify and quantify the stellar deficits, i.e. the flat- 
tened cores, which have long been observed in giant 
galaxies (e.g. King 1978, and references therein), the ra- 
tionale for a physical divide between core-Sersic galax- 
ies and Sersic galaxies also existed long ago. Early-type 
galaxies brighter than Mb = —20.5 mag tend to be 
anisotropic, pressure supported elliptical galaxies with 
boxy isophotes, while the less luminous early- type galax- 
ies tend to have disky isophotes and often contain a ro- 
tating disk of stars (e.g. Carter 1978, 1987; Davies et al. 
1983; Bender et al. 1988; Peletier et al. 1990; Jaffe et 
al. 1994; Faber et al. 1997; Emsellem et al. 2011). In 
addition, and as noted already, the L-a relation changes 
slope at Mb ~ —20.5 mag (Davies et al. 1983; Matkovic 
& Guzman 2005): the core-Sersic galaxies {Mb < —20.5 
mag) define the relation L cx a 5 while the Sersic galax- 
ies (Mb > —20.5 mag) follow the relation L oc a 2 , 
with no discontinuity at Mb = —18 mag. Similarly, 
the elliptical galaxy L-\xq relation, where fiQ is the cen- 
tral surface brightness, also breaks at Mb ~ —20.5 mag 
and remains linear at Mb = —18 mag (e.g. Graham & 
Guzman 2003) 4 . While gas-rich processes and mergers 
are thought to build the Sersic galaxies, relatively gas- 
free galaxy mergers are thought to build the core-Sersic 
galaxies. The coalescence of massive black holes — from 
the pre-merged, gas-poor galaxies — scour out the core of 
the newly merged 'core-Sersic' galaxy (Begelman, Bland- 
ford, & Rees 1980; Ebisuzaki, Makino, & Okumura 1991; 
Merritt, Mikkola, & Szell 2007). 

The above expectation for the 'bent' M^-L relation 
impacts dramatically upon our predictions for black hole 
masses at the centers of other galaxies, has implications 
for the competing formation scenarios proposed to ex- 
plain the coevolution of galactic spheroids and their su- 

3 Sersic galaxies were previously referred to as 'power-law' galax- 
ies before it was realized that their spheroidal component's inner 
light profiles are better described by the curved Sersic model than 
by a power-law (e.g. Trujillo et al. 2004). 

4 As shown, and explained, in Graham 2012c, it is only from 
the use of 'effective' radii and 'effective' surface brightnesses that 
apparent breaks are seen at Mb = — 18 mag in what are actually 
continuous, curved relations. 



permassive black holes, and has ramifications for semi- 
analytic computer codes which try to mimic this coevo- 
lution. Moreover, as noted in Graham (2012a), the real- 
ization of a broken M^-L relation will influence, among 
other things: evolutionary studies of the (black hole)- 
(host spheroid) connection over different cosmic epochs 
(e.g. Cisternas et al. 2011; Li, Ho & Wang 2011; Hiner 
et al. 2012; Zhang, Lu & Yu 2012); predictions for space- 
based gravitational wave detections which use the M^—L 
relation to estimate black hole masses (e.g. Mapelli et al. 
2012, and references therein); and estimates of the black 
hole mass function and mass density based on spheroid 
luminosity functions (e.g. Vika et al. 2009; Li, Wang, & 
Ho 2012). 

Recent studies of the M^-L relation have contin- 
ued to be largely dominated by luminous spheroids with 
Mbh ^ 1O 7 -1O 8 M0, and have continued to fit a straight 
line to the data (e.g. Sani et al. 2011; Vika et al. 2011; 
McConnell et al. 2011a; Beifiori et al. 2012). However, an 
inspection of their M^-L diagrams reveals the onset of 
a steepening relation in their lower luminosity spheroids 
which have Mbh ^ 1O 7 -1O 8 M , as also seen in the dia- 
grams from Graham (2007a) and Giiltekin et al. (2009b). 
In addition, a number of models are actually generating 
'bent' Mbh-Mspheroid relations which show a steepening 
at low masses (e.g. Cirasuolo et al. 2005, their figure 5; 
Khandai et al. 2012, their figure 7; Bonoli, Mayer & Cal- 
legari 2012, their figure 7). Here we build on this (ten- 
tative observational) evidence of a steepening relation at 
low masses by using a sample of 72 galaxies with di- 
rectly measured black hole masses that are tabulated in 
Section 2 and span 10 6 -10 10 M Q . For the first time we 
both identify and quantify the bend in the M^-L rela- 
tion (Section 3). 

Rather than fitting a single quadratic relation to the 
(Mbh, L) distribution of points — which would be in line 
with the use of a log-quadratic relation in the Mhh~n di- 
agram by Graham & Driver (2007a) — we embrace the 
large body of data which shows that Sersic and core- 
Sersic galaxies have distinguishing characteristics sug- 
gestive of a different formation history (e.g. Davies et 
al. 1983; Faber et al. 1997; Graham & Guzman 2003; 
Gavazzi et al. 2005; Ferrarese et al. 2006) and therefore 
fit two Mbh~£ relations, one for each type of galaxy. 

In Section 4 we provide a discussion of some of the 
more salient points arising from the new M^-L relation 
for Sersic and core-Sersic galaxies. In particular, Sec- 
tion 4.1 provides updated black hole mass estimates for 
nearly 50 low-luminosity spheroids, and reveals that they 
occupy the holy grail mass range for intermediate mass 
black holes, ranging from the stellar mass black hole up- 
per limit of ~10 2 M o to the lower-bound of 10 5 -10 6 M Q 
for supermassive black holes. In Section 4.2 we intro- 
duce a significantly modified expression for use in semi- 
analytical models which try to mimic the coevolution 
of supermassive black holes and galaxies. We provide a 
new quadratic cold-gas 'quasar' mode feeding equation 
to match the quadratic (black hole)-to-(host spheroid) 
growth observed in Sersic galaxies built from gas-rich 
processes. Finally, Section 4.3 discusses the observation 
that barred and non-barred disk galaxies currently ap- 
pear to occupy the same region of the M^-L diagram. 

2. DATA 



The M bh -L sphcroid relation 
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We have continued to build on past catalogs of reli- 
able supermassive black hole masses obtained from di- 
rect maser, stellar or gas kinematic measurements (e.g. 
Ferrarese & Ford (2005); Hu (2008); Graham (2008b); 
Graham et al. 2011) by adding 16 galaxies from the lit- 
erature. This gives us an initial sample of 80 galaxies 
with directly measured black hole masses. Since the cat- 
alog in Graham et al. (2011) was prepared, we note that 
several updates over the last two years have resulted in 
the doubling of black hole masses at the high-mass end 
of the Mbh-c diagram (e.g. van den Bosch & de Zeeuw 
2010; Walsh et al. 2010, 2012) and two ten-billion solar 
mass black holes have been reported (e.g. McConnell et 
al. 2012). 

2.1. Galaxy exclusions 

Separate from the above sample of 80 galaxies, we have 
maintained our exclusion of NGC 7457, because its black 
hole mass was derived assuming that this galaxy's excess 
nuclear light was due to an AGN rather than a massive 
nuclear disk and dense star cluster (Balcells et al. 2007; 
Graham & Spitler 2009), and the exclusion of the Sc 
galaxy NGC 2748 due to potential dust complications 
with its estimated black hole mass (Atkinson et al. 2005; 
Hu 2008). 

For the following reasons, our sample of 80 galaxies is 
reduced here by five. The best-fitting parameter set in 
the mass modeling of NGC 5252 by Capetti et al. (2005) 
had resulted in a reduced x 2 value of 16.5, indicative of 
a poor fit to the data and resulting in its exclusion by 
Giiltckin et al. (2009b) . We too now exclude this galaxy, 
which would otherwise appear to have a black hole mass 
that is an order of magnitude too high in the M^-a 
diagram. The inclined, starburst spiral galaxy Circinus 
is also excluded due to its complex kinematics (e.g. For, 
Koribalski & Jarrett 2012, and references therein) and 
location on the other side of the Galactic plane. With a 
Galactic latitude of less than 5 degrees, the suggested B- 
band extinction correction for Circinus is 5 mag or more 
(Schlafly & Finkbeiner 2011) and is considered to be un- 
reliable. The above two galaxies can not be used reliably 
in either the M^-L diagram or the M^-a diagram, and 
they are therefore not listed in our Table [TJ 

Following Giiltekin et al. (2009b), we exclude the 
barred galaxy NGC 3079. Although NGC 3079 has a well 
determined black hole mass (Kondratko et al. 2005), the 
observed stellar velocity dispersion drops rapidly from a 
central value of ^150 km s _1 to 60 km s _1 or less within 
just 1 kpc (Shaw, Wilkinson, & Carter 1993). The dy- 
namics in this Sc galaxy's boxy peanut-shaped bulge are 
known to be affected by bar streaming motions (Merri- 
field & Kuijken (1999; Veilleux, Bland-Hawthorn, k Ce- 
cil 1999) which were not modeled by Shaw et al. (1993). 
Such additional motions are a general concern in barred 
galaxies, as detailed in Graham et al. (2011), and it is 
not clear what velocity dispersion should be used for the 
bulge component of this galaxy 5 . The bulge luminosity 
of the Milky Way is uncertain due to dust in our Galaxy's 
disc, and is therefore not included in our Mbh-£ S pheroid 
diagram. Finally, we also exclude M32 which belongs to 

5 While this galaxy could be used in the Mt,h-L S pheroid diagram, 
we do not so as to be consistent with our galaxy sample in the final 
Mbh-f diagram. 



the rare 'compact elliptical' galaxy class and is thus not 
representative of the majority of galaxies. Such galaxies 
are thought to be heavily stripped of an unknown frac- 
tion of their stars (e.g. Bekki et al. 2001), and as such 
this galaxy may bias our analysis. We do however note 
that M32, and the Milky Way, have been included in 
a preliminary M^-a diagram (Figure [T|). We therefore 
mention here that we have used a host bulge velocity dis- 
persion of 55 km s _1 from M32's bulge stars outside of 
M32's core region; this value is lower than the commonly 
used central aperture values of 72-75 km s _1 which are 
biased high by the dynamics around the central black 
hole (Igor Chilingarian, priv. comm.). This leaves us, 
thus far, with a sample of (80 — 5 —) 75 galaxies. For 
reasons discussed at the start of Section |3l an additional 
three galaxies (NGC 1316, NGC 3842 and NGC 4889) 
are excluded from some of the analyses. 

2.2. Galaxy distances and magnitudes 

For the bulk of the sample, we have used the dis- 
tance moduli from Tonry et al. (2001) after first de- 
creasing these values by 0.06 mag and thereby reducing 
both the galaxy distances and the black hole masses by 
~3%. This small adjustment arises from Blakeslee et 
al.'s (2002, their Section 4.6) recalibration of the sur- 
face brightness fluctuation method based on the final 
Cepheid distances given by Freedman et al. (2001, with 
the metallicity correction). In addition to listing each 
galaxy's distance and velocity dispersion, primarily taken 
from HyperLeda 6 (Paturel et al. 2003), Table [T] presents 
their observed (i.e. uncorrected) B-band magnitude, Bt, 
from the Third Reference Catalogue of Bright Galaxies 
(de Vaucouleurs et al. 1991, hereafter RC3) and their to- 
tal -ftT s -band magnitude as provided by the Two Micron 
All Sky Survey (2MASS) 7 (Jarrett et al. 2000). Due to 
the bright star near the center of NGC 2974, we used the 
K-b&nd magnitude from Cappellari et al. (2006) rather 
than the 2MASS value for this galaxy. 

Both the observed B- and if s -band magnitudes need 
to be corrected for Galactic extinction, which we have 
taken from Schlegel et al. (1998), as provided by the 
NASA/IPAC Extragalactic Database (NED) 8 and in- 
cluded in Table 1 as As and Ak for ease of reference. 
We additionally provide K-corrections for the 2MASS 
-Kg-band magnitudes, derived from each galaxy's helio- 
centric redshift z (taken from NED) and total J — K 
color (taken from 2MASS) coupled with the prescription 
from Chilingarian, Melchior & Zolotukhin (2010) which 
is available online 9 . In passing, we note that for the low 
redshifts associated with our galaxy sample, this 
correction is approximately —2.12;, — 2.2z, and — 2.5z for 
the elliptical, SO and Sa, and later spiral galaxy types, 
respectively. Without reliable B — R or B — I colors, the 
£?-band K-corrections tabulated in Table 1 are such that 
Kb equals 4z, 3.5z, 3z, and 2.5z for the elliptical, SO, 
Sa, and later spiral galaxy types, respectively. As can 
be seen in Table [TJ these corrections are minor and they 
have had no significant impact on the final results. Fi- 
nally, in our analysis we corrected all galaxy magnitudes 

6 http://leda.univ-lyonl.fr 

7 www.ipac.caltech.edu/2mass 

8 http://nedwww.ipac.caltech.edu 

9 http://kcor.sai.msu.ru 
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for cosmological redshift dimming. 



and 



2.3. Bulge magnitudes 

In addition to the above three corrections that were 
applied in the analysis rather than in Table [TJ the tabu- 
lated disk galaxy magnitudes were adjusted for two other 
factors: their observed, total magnitudes were converted 
into dust- corrected, bulge magnitudes. Rather than do 
this galaxy by galaxy (e.g. Grootes et al. 2012), which 
would require careful bulge/disk decompositions, we can 
take advantage of our large sample size and employ a 
mean statistical correction. While this will result in in- 
dividual bulge magnitudes not being exactly correct, the 
ensemble average correction will be correct. Given that 
dust has an order of magnitude less impact in the K s - 
band than in the £?-band, we can expect there to be less 
scatter in our AVband Mbh~L relation, at least for the 
Sersic sample which contains the bulk of the bulges. 

Using the relation M = — 2.51og(L), the observed and 
the intrinsic (i.e. dust-corrected) bulge and disk lumi- 
nosities are related as follows. 

£bulge,obs/£bulge,int = 10~ AMbul s°/ 2 - J , (1) 
-^disk.obs/idiskant = 10~ AMdlsk / 2 - 5 , (2) 

r, , , _|_ T.j. i , — 1f)- M gai«y,obs/2.5 /Q\ 
-^bulgc,obs T -^disk.obs — J-U ; V° I 

-£'bulge,int/-£'disk,int = B/D, (4) 

where AMbuigc and AMdisk are the differences between 
the observed and intrinsic magnitudes (given below). 
The above equations can be combined to give the use- 
ful expression 



Mi 



bulge, int " 



- M, 



galaxy, obs 



-2.5 log 



10 -AA/ bulg c/2.5 _|_ 



^Q-AM diBk /2.5 



(B/D) 



5) 



The average dust-corrected bulge-to-disk (B/D) flux 
ratio is well known to be a function of both galaxy mor- 
phological type and passband, varying mainly due to the 
bulge luminosity (e.g. Yoshizawa & Wakamatsu 1975; 
Trujillo et al. 2001). Graham & Worley (2008) presented 
over 400 if- band, B/D flux ratios as a function of disk 
galaxy morphological type. Given that extinction due 
to dust is still an issue at near-infrared wavelengths, al- 
beit notably less in the if-band than the 7i-band, these 
flux ratios were corrected by Graham & Worley (2008) 
for dust following the prescription given by Driver et al. 
(2008). Most studies have not made this important cor- 
rection which effectively increases the observed B/D flux 
ratio after accounting for the light from the far-side of 
the bulge which is obscured by the centrally-concentrated 
dust in the disk. Here we have slightly adjusted Graham 
& Worley's (2008) B/D values by combining the SO and 
SO/a classes as one in our Tableland ensuring that they 
do not have a smaller B/D ratio than the Sa galaxies. 
By combining these ratios with the average, inclination- 
dependent AMbuigc and AMdisk values below, one can 
derive the expected 'intrinsic' bulge magnitude from the 
observed (dust-dimmed) galaxy magnitude and the incli- 
nation of the disk. Specifically, for AMb u i ge and AMdisk 
we have used the expressions from Driver et al. (2008) 
such that 



AM, 



bulge 



= M 



bulge, obs" 



-M, 



bulge, intrin 



= b 1 -\-b 2 [1 — cos(z)] fc 



(6) 



AMdisk = M&sc.obs - Mdisc.intrin = dl + d 2 [l - COs(z)] 3 , 

(7) 

where the coefficients in the above equations depend 
on the passband used and are provided in Driver et al. 
(2008). The cosine of the disk inclination angle i, such 
that i = degrees for a face-on disk, is roughly equal 
to the minor-to-major axis ratio, b/a, at large radii. For 
our galaxy sample, these values have been provided in 
the final column of Table [1] 

We have used a conservative (small) uncertainty of 5% 
for the velocity dispersions, and assigned a typical un- 
certainty of 0.25 mag to the elliptical galaxy magnitudes. 
Regarding the disk galaxies, the observed range of bulge- 
to-disk flux ratios (Graham & Worley 2008) has a la 
scatter equal to a factor of ~2 for any given disk galaxy 
morphological type. We therefore assign a notably larger 
uncertainty of 0.75 mag to our bulge magnitudes. While 
this may at first sound worryingly high, we again note 
that the sample average shift from a disk galaxy magni- 
tude to a bulge magnitude will be more accurate, thereby 
much less affecting the recovery of the M^-L relation. 
It is only now that the disk galaxy sample size N = 44 
is sufficiently large (coupled with 9 elliptical galaxies in 
the Sersic sample) that we can use this methodology be- 
cause the uncertainty in the disk-galaxy sample average 
correction scales with 1/y/N. 

2.4. Galaxy core type 

As described in the Introduction, several global prop- 
erties and scaling relations of galaxies depend on their 
core type. We have therefore identified galaxies with or 
without a partially depleted stellar core. Such cores, with 
shallow light profiles, are thought to have formed during 
'dry' galaxy merger events in which a binary supermas- 
sive black hole's orbit decays by gravitationally ejecting 
stars from the center of the newly merged galaxy (e.g. 
Begelman et al. 1980; Ebisuzaki et al. 1991; Graham 
2004; Merritt et al. 2007). However Sersic galaxies that 
have n < 2-3 may also have a resolved inner, negative 
logarithmic light-profile slope (for the spheroid) which 
is less steep than 0.3, although they have no depleted 
core relative to their outer profile. Graham et al. (2003) 
warned that some of these galaxies may be misidentificd 
as 'core' galaxies and Dullo & Graham (2012) have shown 
that this has occurred ~20% of the time in the past. We 
therefore have a preference for cores identified as a deficit 
relative to the inward extrapolation of the spheroid's 
outer Sersic profile (e.g. Trujillo et al. 2004; Ferrarese 
et al. 2006; Dullo & Graham 2012). Such galaxies are 
referred to as 'core-Sersic' galaxies, while those with no 
deficit or instead additional nuclear components are re- 
ferred to as 'Sersic' galaxies (e.g. Graham et al. 2003; 
Graham & Guzman 2003; Balcells et al. 2003; Trujillo et 
al. 2004). 

For most of the luminous galaxies we have been able 
to identify from suitably high-resolution images whether 
or not they possess a partially-depleted core, although 
some are at too great a distance to know. At the low- 
luminosity end of the relation, dusty nuclei can make this 
task more challenging. When no core designation was 
available or possible from the literature (see column 4 
of Table [TJ we used the velocity dispersion to help us 
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TABLE 1 
Black hole / Galaxy data. 



Galaxy 


Type 


Dist core a A^bh 


B'j' 


K s 


A s A K K s K Ks 


R-25 


M s 


M K , 






[Mpc] km s- 1 [1O 8 M ] 


[mag] 


[mag] 


[mag] [mag] [mag] [mag] 




[mag] 


[mag] 


(1) 


(2) 


(3) (4) (5) (6) 


(7) 


(8) 


(9) (10) (11) (12) 


(13) 


(14) 


(15) 



A1836, BCG El 


158.0 [3a] 


y? 


309 


[5a] 


A3565, BCG El 


40.7 [3a] 


y? 


335 


[5b] 


IC 1459 


E3 


28.4 


y? 


306 




NGC 821 


E6 


23.4 


n [4a] 


200 




NGC 1399 


El 


19.4 


y [4b] 


329 




NGC 2974 


E4 


20.9 


n [4c] 


227 




NGC 3377 


E5 


10.9 


n [4d] 


139 




NGC 3379 


El 


10.3 


y [4b] 


209 




NGC 3608 


E2 


22.3 


y [4b] 


192 




NGC 3842 


E3 


98.4 [6h] 


y [4b] 


270 


[6h] 


NGC 4261 


E2 


30.8 


y [4c] 


309 




NGC 4291 


E3 


25.5 


y [4b] 


285 




IN 4o / 4 




17.9 


y [4f] 


296 


[6j] 


NGC 4473 


E5 


15.3 


n [4b] 


179 




NGC 4486 


El 


15.6 


y [4f] 


334 




NGC 4486a 


E2 


17.0 [3b] 


n [4f] 


110 


[5c] 


NGC 4552 


E0 


14.9 


y [4b] 


252 




NGC 4621 


E5 


17.8 


n [4f] 


225 




NGC 4649 


E2 


16.4 


y [4b] 


335 




NGC 4697 


E6 


11.4 


n [4g] 


171 




NGC 4889 


E4 


103.2 [6h] 


y [4g] 


347 


[6h] 


NGC 5077 


E3 


41.2 [3a] 


y [4h] 


255 




NGC 5576 


E3 


24.8 


n [4b] 


171 




NGC 5813 


El 


31.3 


y [4b] 


239 




NGC 5845 


E4 


25.2 


n [4h] 


238 




NGC 5846 


E0 


24.2 


y [4i] 


237 




NGC 6086 


E3 


138.0 [3a] 


y [4j] 


318 


[6p] 


NGC 6251 


E2 


104.6 [3a] 


y? [6q] 


311 




NGC 7052 


E4 


66.4 [3a] 


y [4k] 


277 




NGC 7768 


E2 


112.8 [6h] 


y [41] 


257 


[6h] 


Cygnus A 


Sa? 


232.0 [3a] 


y? 


270 




IC 2560 


SBb 


40.7 [3a] 


n? [4m] 


144 


[5d] 


NGC 224 


Sb 


0.74 


n [4n] 


170 




NGC 253 


SBc 


3.5 [3c] 


n [4o] 


109 


[5c] 


NGC 524 


SO 


23.3 


y [4c] 


253 




NGC 1023 


SB0 


11.1 


n [4g] 


204 




NGC 1068 


SBb 


15.2 [3a] 


n [4p] 


165 


[5f] 


NGC 1194 


SO 


53.9 [3a] 


n? 


148 


[5g] 


NGC 1300 


SBbc 


20.7 [3a] 


n [4q] 


229 




NGC 1316 


SB0 


18.6 [3d] 


? 


226 




NGC 1332 


SO 


22.3 


y? 


320 




NGC 2273 


SBa 


28.5 [3a] 


n [4r] 


145 


[5g] 


NGC 2549 


SB0 [2a] 


12.3 


n [4q] 


144 




NGC 2778 


SB0 [2b] 22.3 


n [4d] 


162 




NGC 2787 


SB0 


7.3 


n [4s] 


210 




NGC 2960 


Sa? 


81.0 [3c] 


n? 


166 


[5g] 


NGC 3031 


Sab 


3.8 


n [4t] 


162 




NGC 3115 


SO 


9.4 


n [4a] 


252 




NGC 3227 


SBa 


20.3 [3a] 


n [4c] 


133 




NGC 3245 


SO 


20.3 


n [4c] 


210 




NGC 3368 


SBab 


10.1 


n [4u] 


128 




NGC 3384 


SB0 


11.3 


n [4a] 


148 




NGC 3393 


SBa 


55.2 [3a] 


n [4m] 


197 




NGC 3414 


SO 


24.5 


n [4c] 


237 





[6a] 



39±| 

[6a] 
24ii° [6b] 
0.39i«;26 [6c] 
4.7i°'6 [6d] 
1.7i°; 2 2 [6c] 
0.77tg;{£ [6f] 
4.0±i;° [6g] 

2-oiJi 

97±f 6 [6h] 
5.0±i"g [6i] 

" [6c] 
[6j] 



30 Ellipticals 



3.3i°§ 
9.0i°j 

i-2±g:S [6c] 

58±| | [6k] 
0.13+°-°! [61] 
4.7i°; 5 5 [6c] 
3.9i°-j [6c] 
47iJ° [6m] 
" [6c] 
210il|g [6h] 
7A±t 7 [6n] 



1 8+ - 2 
1 - 8 -0.1 



1- e±8:2 tH 

6.8i°;7 [6e] 

2- 6i° j [6c] 
llii [6c] 
37il? [6p] 
5.9i|° [6q] 

3- 7i?;l [6r] 
13±* [6h] 



14.56 
11.61 
11.61 
11.67 
10.55 
11.87 
11.24 
10.24 
11.70 
12.78 
11.41 
12.43 
10.09 
11.16 
09.59 
13.20 [7a] 
10.73 
10.57 
9.81 
10.14 
12.53 
12.38 
11.85 
11.45 
13.50 
11.05 
13.79 
13.64 [7b] 
13.40 [7b] 
13.24 



25± 7 7 



45 Bulges 



. [6s] 
0.044i°;°H [6t] 
1.4tg;» [6u] 
o.ioi°;J° [6v] 

0.42; 



i 0.04 



-0.04 
0.084i° 0°| [6y] 

0-66i°;| [6z] 

; 35 [6aa] 
i k+0.75 
± -°-0.8 
14.5i| [6ac 
0. 

,004 

14+ 002 
0.15T°° 9 [6c] 
0.40i° « [ 6 ad] 



0.731 

[6ab] 

+2 

0.083+°';°°! [6 Z 
[6w] 



0,05 
0.005 [ 6z 



q 117 +b.oos 



[6af] 



0.74+°;£ [6ae] 

o O + 10.0 

0-14^06 [6ag] 
2.0±g | [6ah] 
0.073±°;°i| [6af 
0-i7±g;8i [6c] 
0.34±°;{g [6aj] 
2.4±g:l [6e] 



17.04 

12.53 [7b] 

4.36 

8.04 

11.30 

10.35 

9.61 

13.83 

11.11 

9.42 

11.25 

12.55 

12.19 

13.35 

11.82 

13.29 [7b] 

7.89 

9.87 

11.10 

11.70 

10.11 

10.85 

13.09 [7b] 

11.96 



9.993 

7.502 

7.502 

7.900 

6.306 

8.00 [8a] 

7.441 

6.270 

8.096 

9.082 

7.263 

8.417 

6.222 

7.157 

5.812 

9.012 

6.728 

6.746 

5.739 

6.367 

8.407 

8.216 

7.827 

7.413 

9.112 

6.949 

9.973 

9.026 

8.574 

9.335 

10.276 

8.694 

0.984 

3.772 

7.163 

6.238 

5.788 

9.758 

7.564 

5.587 

7.052 

8.480 

8.046 

9.514 

7.263 

9.783 

3.831 

5.883 

7.639 

7.862 

6.320 

6.750 

9.059 

7.981 



0.277 
0.265 
0.265 
0.474 
0.056 
0.235 
0.147 
0.105 
0.090 
0.093 
0.078 
0.158 
0.174 
0.123 
0.096 
0.102 
0.177 
0.143 
0.114 
0.131 
0.041 
0.210 
0.136 
0.246 
0.230 
0.237 
0.162 
0.377 
0.522 
0.167 

1.644 
0.410 
0.268 
0.081 
0.356 
0.262 
0.145 
0.330 
0.130 
0.090 
0.141 
0.305 
0.282 
0.090 
0.565 
0.193 
0.346 
0.205 
0.098 
0.108 
0.109 
0.115 
0.325 
0.106 



0.024 
0.023 
0.023 
0.040 
0.005 
0.020 
0.013 
0.009 
0.008 
0.008 
0.007 
0.013 
0.015 
0.010 
0.008 
0.009 
0.015 
0.012 
0.010 
0.011 
0.004 
0.018 
0.012 
0.021 
0.020 
0.020 
0.014 
0.032 
0.044 
0.014 

0.140 
0.035 
0.023 
0.007 
0.030 
0.022 
0.012 
0.028 
0.011 
0.008 
0.012 
0.026 
0.024 
0.008 
0.048 
0.016 
0.029 
0.017 
0.008 
0.009 
0.009 
0.010 
0.028 
0.009 



0.15 
0.05 
0.02 
0.02 
0.02 
0.03 
0.01 
0.01 
0.02 
0.08 
0.03 
0.02 
0.01 
0.03 
0.02 
0.00 
0.00 
0.01 
0.01 
0.02 
0.09 
0.04 
0.02 
0.03 
0.02 
0.02 
0.13 
0.10 
0.06 
0.11 

0.17 
0.02 
0.00 
0.00 
0.03 
0.01 
0.01 
0.05 
0.01 
0.02 
0.02 
0.02 
0.01 
0.02 
0.01 
0.05 
0.00 
0.01 
0.01 
0.02 
0.01 
0.01 
0.04 
0.02 



-0.07 
-0.03 
-0.01 
-0.01 
-0.01 
-0.01 
-0.01 
-0.01 
-0.01 
-0.04 
-0.02 
-0.01 
-0.01 
-0.02 
-0.01 
-0.00 
-0.00 
-0.00 
-0.01 
-0.01 
-0.04 
-0.02 
-0.01 
-0.01 
-0.01 
-0.01 
-0.05 
-0.04 
-0.03 
-0.05 



-0.07 0.13 
-0.02 0.20 
-0.00 0.49 
-0.00 0.61 
-0.02 0.01 
-0.00 0.47 
-0.01 0.07 
-0.03 0.25 
-0.01 0.18 
-0.02 0.15 
-0.01 0.51 
-0.01 0.12 
-0.01 0.48 
-0.02 0.14 
-0.00 0.19 
-0.03 0.16 
-0.00 0.28 
-0.01 0.47 
-0.01 0.17 
-0.01 0.26 
-0.01 0.16 
-0.01 0.34 
-0.02 0.04 
-0.01 0.14 



-21.72 
-21.71 
-21.37 
-20.65 
-20.95 
-19.97 
-19.09 
-19.93 
-20.13 
-22.28 
-21.11 
-19.76 
-21.35 
-19.89 
-21.47 
-18.05 
-20.31 
-20.83 
-21.38 
-20.28 
-22.59 
-20.91 
-20.26 
-21.28 
-18.74 
-21.11 
-22.08 
-21.84 
-21.24 
-22.20 

-20.80 
-19.02 
-18.67 
-17.33 
-20.15 
-19.82 
-19.45 
-19.58 
-18.05 
-21.33 
-20.37 
-19.33 
-18.24 - 
-17.79 
-17.41 
-20.79 
-19.59 
-19.88 
-19.87 
-19.37 
-19.14 
-19.05 
-20.21 
-19.40 



-26.25 
-25.65 
-25.50 
-24.02 
-25.17 
-23.66 
-22.78 
-23.83 
-23.68 
-26.02 
-25.24 
-23.66 
-25.08 
-23.83 
-25.19 
-22.15 
-24.16 
-24.52 
-25.37 
-23.96 
-26.80 
-24.94 
-24.19 
-25.12 
-22.95 
-25.02 
-25.93 
-26.25 
-25.68 
-26.11 

-25.81 
-22.78 
-21.82 
-21.42 
-23.61 
-23.01 
-23.47 
-22.91 
-21.91 
-24.68 
-23.74 
-22.66 
-21.44 
-21.15 
-20.98 
-23.69 
-22.56 
-23.00 
-22.75 
-22.61 
-22.17 
-22.47 
-23.56 
-22.87 
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TABLE 1 
cord. 





lypc 


Dist 


core 




lvl bh 





K 

"■s 


A r, 


Ax 


K s 


K „ 
^K s 


^25 


A/T ^ 
1V1.B 


A/T „ 






[MdcI 




km s ^ 


[10 s Mml 


fmad 


[masl 


[ma£rl 


[mag] [mag] 


[masl 




[mael 


[matd 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


NGC 3489 


SBO 


11.7 


n [4q] 


105 


n nrro + 0.008 xa -i 
°- 058 -0.008 [ 6ai l 


11.12 


7.370 


0.072 


0.006 


0.01 


—0.01 


0.24 


— 18.69 


—21.88 


NGC 3585 


SO 


19.5 


n [4q] 


206 


o i +1.4 i 

3.il .6 M 


10.88 


6.703 


0.276 


0.023 


0.02 


—0.01 


0.26 


—20.27 


-23.69 


NGC 3607 


SO 


22.2 


n [4e] 


224 


l-3Io.5 M 


10.82 


6.994 


0.090 


0.008 


0.01 


—0.01 


0.30 


—20.47 


-23.68 


NGC 3998 


SO 


13.7 


y? [4v] 


305 


1 +2 . X r- 1 1 

8.1_! 9 |6ak] 


11.61 


7.365 


0.069 


0.006 


0.01 


—0.01 


0.08 


— 18.41 


—22.21 


NGC 4026 


SO 


13.2 


n [4q] 


178 


1 0+O.6 Xn l 

1.8I 3 [60] 


11.67 


7.584 


0.095 


0.008 


0.01 


—0.01 


0.61 


— 18.88 


—22.11 


NGC 4151 


SBab 


20.0 [3aJ 


r a l 
n [4eJ 


156 


0.651 07 [bal] 


11.50 


7.381 


0.119 


0.010 


0.01 


—0.01 


0.15 


— 19.24 


—22.60 


NGC 4258 


SBbc 


7.2 [3f] 


r a l 
n [4wJ 


134 


0.39_ 01 [bam] 


9.10 


5.464 


0.069 


0.006 


0.00 


—0.00 


0.41 


— 17.95 


—21.76 


NGC 4342 


SO 


23.0 [3g] 


r a l 
n [4xJ 


253 


4.5_j g ]6an] 


13.41 


9.023 


0.088 


0.008 


0.01 


—0.01 


0.33 


— 17.99 


—21.73 


NGC 4388 


Sb 


17.0 [3b] 


o r a i 
n? [4y] 


107 [5g] 


°- 075 -0,OO2 M 


11.76 


8.004 


0.143 


0.012 


0.02 


—0.02 


0.64 


—20.07 


—23.65 


NGC 4459 


SO 


15.7 


n [4f] 


178 


n f n + . 1 3 X f ll 

0.68_y ^ [6ad] 


11.32 


7.152 


0.199 


0.017 


0.01 


—0.01 


0.12 


-19.15 


—22.73 


NGC 4564 


SO [2c] 


14.6 


n [4f] 


157 


0.60±° : {g [6c] 


12.05 


7.937 


0.151 


0.013 


0.01 


-0.01 


0.38 


-18.49 


-21.87 


NGC 4594 


Sa 


9.5 


y [6ao] 


297 [6ao] 


6.4l«;4 [6ao] 


8.98 


4.962 


0.221 


0.019 


0.01 


-0.01 


0.39 


-20.71 


-23.86 


NGC 4596 


SBO 


17.0 [3b] 


n [4z] 


149 


°- 79+ 33 [ 6ad l 


11.35 


7.463 


0.096 


0.008 


0.02 


-0.01 


0.13 


-19.20 


-22.60 


NGC 4736 


Sab 


4.4 [3i] 


n? 


104 


0.060±g"git [6ap] 


8.99 


5.106 


0.076 


0.007 


0.00 


-0.00 


0.09 


-18.37 


-21.55 


NGC 4826 


Sab 


7.3 


n? 


91 


0-016t° °°t [6ap] 


9.36 


5.330 


0.178 


0.015 


0.00 


-0.00 


0.27 


-19.36 


-22.47 


NGC 4945 


SBcd 


3.8 [3i] 


n? [4aa] 


100 


0-014j:°-^ [6aq] 


9.30 


4.483 


0.762 


0.065 


0.00 


-0.00 


0.72 


-16.39 


-20.60 


NGC 5128 


SO 


3.8 [3i] 


n? [4ab] 


120 


0.45±g"i5 [6ar] 


7.84 


3.942 


0.496 


0.042 


0.01 


-0.00 


0.11 


-19.84 


-22.87 


NGC 6264 


Sa? 


146.3 [3a] 


n? 


159 [5g] 


0.3051°-°°* [6z] 


15.42 [7b] 11.407 0.280 


0.024 


0.10 


-0.07 


0.15 


-20.04 


-23.47 


NGC 6323 


SBab [2d] 112.4 [3a] 


n? 


159 [5g] 


o-io±g:Sgi [ 6Z ] 


14.85 [7b] 10.530 0.074 


0.006 


0.06 


-0.05 


0.48 


-20.01 


-23.44 


NGC 7582 


SBab 


22.0 [3a] 


n [4ac] 


156 


0.55l°;?| [6as] 


11.37 


7.316 


0.061 


0.005 


0.01 


-0.01 


0.38 


-19.77 


-22.94 


UGC 3789 


SBab 


48.4 [3a] 


n? 


107 [5g] 


0.1081°;°°* [6z] 


13.30 [7b] 9.510 


0.280 


0.024 


0.03 


-0.04 


0.06 


-19.48 


-22.46 












3 extra 




















NGC 3079 


SBc 


20.7 [3a] 


n? [4ad] 


60-150 [5h] 0.024i2;gi2 [6at] 


11.54 


7.262 


0.049 


0.004 


0.01 


-0.00 


0.74 


-17.80 


-21.86 


Milky Way SBbc 


0.008 [6au] 


n [4ae] 


100 [5i] 


0.043l° °° 4 [6au] 




















M32 


E? [2e] 


0.79 


n [4ae] 


55 [5j] 


0.024l°;°° 5 5 [6av] 


9.03 


5.095 


0.268 


0.023 


0.00 


-0.00 


0.13 


-15.73 


-19.42 



Note. — Column 1: Galaxy name. Column 2: Basic morphological type, primarily from NED". Column 3: Distance, primarily from 
Tonry et al. (2001) and corrected according to Blakcslcc ct al. (2002). Column 4: Presence of a partially depleted core. The addition of a 
question mark is used to indicate that this classification has come from the velocity dispersion criteria mentioned in the text. Column 5: 
Velocity dispersion primarily from HypcrLeda 6 (Paturel et al. 2003). Column 6: Black hole mass, adjusted to the distances in column 3. 
Column 7: Bt is the total (observed) magnitude in the B system from the RC3, unless otherwise noted. Column 8: K s is the total 
2MASS A" 3 -band magnitude. Columns 9 & 10: The B and i-T-band Galactic extinction from Schlegcl ct al. (1998). Columns 11 & 12: 
The B and iiT-band K-correction. Column 13: R25, taken from the RC3, is the (base 10) logarithm of the major-to-minor (a/6) isophotal 
axis ratio at /is =25 mag arcsec -2 . Columns 14 & 15: Corrected, absolute, B- and if-band spheroid magnitude. References: 2a = 
Krajnovic ct al. (2009); 2b = Rix et al. (1999); 2c = Trujillo ct al. (2004); 2d = Jarrett et al. (2000) 2MASS image; 2e = Graham 
(2002); 3a = NED (Virgo + GA + Shapley)-corrected Hubble flow distance; 3b = Jerjen et al. (2004); 3c = Rekola et al. (2005); 3d 
= Madore et al. (1999); 3e = Violette Impellizzeri et al. (2012); 3f = Herrnstein et al. (1999); 3g = Mei et al. (2007), Blakeslee et al. 
(2009), Blom ct al. (2012), Bogdan et al. (2012b); 3h = Herrmann ct al. (2008); 3i = Karachentsev et al. (2007); 4a = Ravindranath et 
al. (2001); 4b = Dullo & Graham (2012); 4c = Lauer et al. (2005); 4d = Rest et al. (2001); 4e = Richings et al. (2011); 4f = Ferrarese 
et al. (2006); 4g = Faber et al. (1997); 4h = Trujillo et al. (2004); 4i = Forbes, Brodie, & Huchra (1997); 4j = Laine et al. (2003); 4k = 
Quillcn et al. (2000); 41 = Grillmair et al. (1994); 4m = Munoz Marin et al. (2007); 4n = Corbin, O'Neil, & Rickc (2001); 4o = Kornei 
& McCrady (2009); 4p = Macchetto et al. (1994); 4q = Graham (2012b); 4r = Malkan, Gorjian & Tarn (1998); 4s = Erwin ct al. (2003); 
4t = Fisher & Drory (2008); 4u = Nowak et al. (2010); 4v = Gonzalez Delgado et al. (2008); 4w = Pastorini et al. (2007); 4x = van den 
Bosch, Jaffe, & van der Marel (1998); 4y = Pogge & Martini (2002); 4z = Gerssen, Kuijken, & Merrifield (2002); 4aa = Marconi ct al. 
(2000); 4ab = Radomski et al. (2008); 4ac = Wold & Galliano (2006); 4ad = Cecil et al. (2001); 4ae = Graham & Spitler (2009); 5a = 
Dalla Bonta et al. (2009); 5b = Smith et al. (2000); 5c = De Francesco et al. (2006); 5d = Cid Fernandes et al. (2004); 5e = Oliva et 
al. (1995); 5f = Nelson & Whittle (1995); 5g = Greene et al. (2010); 5h = Shaw, Wilkinson, & Carter (1993); 5i = Merritt & Ferrarese 
(2001a); 5j = Lucey et al. (1997), Chilingarian (2012, in prep.); 6a = Dalla Bonta et al. (2009); 6b = Cappellari et al. (2002), stellar 
dynamical measurement; 6c = Gebhardt et al. (2003), Giiltckin ct al. (2009b); 6d = Houghton ct al. (2006), Gebhardt et al. (2007); 6c 
= Preliminary values determined by Hu (2008) from Conf. Proc. figures of Cappellari et al. (2008); 6f = Copin et al. (2004); 6g = van 
den Bosch & de Zeeuw (2010); 6h = McConncll ct al. (2012); 6i = Ferrarese et al. (1996); 6j = Walsh et al. (2010); 6k = Gebhardt ct al. 
(2011); 61 = Nowak et al. (2007); 6m = Shen & Gebhardt (2009); 6n = De Francesco ct al. (2008); 60 = Giiltckin ct al. (2009a); 6p = 
McConnell et al. (2011b); 6q = Ferrarese & Ford (1999); 6r = van der Marel & van den Bosch (1998); 6s = Tadhunter et al. (2003); 6t 
= Ishihara et al. (2001), Nakai et al. (1998); 6u = Bacon et al. (2001), Bender et al. (2005); 6v = Rodriguez-Rico et al. (2006), a factor 
of 2 uncertainty has been assigned here; 6w = Krajnovic ct al. (2009); 6x = Bower ct al. (2001); 6y = Lodato & Bertin (2003); 6z = 
Kuo et al. (2011); 6aa = Atkinson ct al. (2005); 6ab = Nowak et al. (2008); 6ac = Rush et al. (2011); 6ad = Sarzi et al. (2001); 6ae = 
Devereux et al. (2003); 6af = Emsellem et al. (1999); 6ag = Davies ct al. (2006), Hicks & Malkan (2008); 6ah = Barth et al. (2001); 6ai 
= Nowak ct al. (2010); 6aj = Kondratko et al. (2008); 6ak = Walsh ct al. (2012); 6al = Onken et al. (2007), Hicks & Malkan (2008); 
6am = Herrnstein ct al. (1999); 6an = Crctton & van den Bosch (1999), Valluri ct al. (2004); 6ao = Jardcl ct al. (2011), NGC 4594 is 
a core-Sersic galaxy with an AGN; 6ap = Kormendy et al. (2011), priv. value from K. Gebhardt; 6aq = Grccnhill et al. (1997); 6ar = 
Neumayer (2010); 6as = Wold et al. (2006); 6at = Trotter et al. (1998), Yamauchi et al. (2004), Kondratko et al. (2005); 6au = Gillessen 
et al. (2009); 6av = Verolme et al. (2002); 7a = B T (Johnson), Gavazzi et al. (2005); 7b = The RC3's m B value (which they reduced 
to the B-r system); 8a = Cappellari et al. (2006). 

a http: / / ncdwww.ipac.caltech.edu 
^http^/leda .univ-lyonl.fr 
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TABLE 2 

Average dust-corrected (intrinsic) bulge-to-disk flux 

RATIOS 



Galaxy <log(B/D)> <log(B/D)> 
Type B-band K-band 



SO/SOa 


-0.29 


-0.31 


Sa 


-0.29 


-0.34 


Sab 


-0.39 


-0.54 


Sb 


-0.87 


-0.60 


Sbc 


-1.13 


-0.82 


Sc 


-1.28 


-1.06 


Scd 


-1.55 


-1.23 



Adapted from Graham & Worley (2008). 



assign the core type. This meant that for seven galax- 
ies with a > 270 km s -1 , all but one of which actu- 
ally have a > 300 km s _1 , we classified them as core- 
Sersic galaxies. For 12 galaxies with a < 166 km s _1 wc 
have classified them as Sersic galaxies. This approach 
relates to the trend in which low-luminosity spheroids, 
often found in disk galaxies, do not possess partially de- 
pleted cores while luminous spheroids do (e.g. Nieto et 
al. 1991; Ferrarese et al. 1994; Faber et al. 1997). For 
the remaining galaxies, their core type was determined 
from high- resolution images (see column 4 of Tabic [1]). 

While Faber et al. (1997) reported that the peculiar 
galaxy NGC 1316 (Fornax A) has a partially-depleted 
core, we feel that a proper bulge/disk/bar decomposi- 
tion — which is beyond the scope of this paper — might 
be required for this barred SO galaxy. This may reveal 
a bulge with a small Sersic index and a shallow inner 
profile slope, rather than an actual depleted core. Al- 
though NGC 1316 is a merger remnant (Deshmukh et al. 
2012, and references therein) — which may support the 
presence of a depleted core if it had been a dry merger 
event — the presence of much dust and gas, plus the ex- 
istence of a bar, is indicative of a minor, wet merger. We 
note that no other barred galaxy in our sample has a core 
which is partially depleted of stars, and to be safe we felt 
it best to exclude this galaxy from our final analysis until 
a more detailed investigation of its stellar distribution is 
performed. 

3. ANALYSIS AND RESULTS 

Our initial investigation of the M^-a diagram uses 
all 75 galaxies (including NGC 1316) plus M32 and the 
Milky Way's bulge to give a total sample of 77 galaxies. 
It is apparent that there is something unusual with the 
two ultramassive black holes reported for NGC 3842 and 
NGC 4889 by McConnell et al. (2011a). Figure □ reveals 
that they may be are outliers biasing the M^-a relation 
defined by the remaining data. Performing a symmetri- 
cal linear regression on the core-Sersic galaxies, including 
NGC 3842 and NGC 4889, yields a slope of roughly 7.0±1 
(see Table |3|) . Uncertain as to whether there is an error 
in the black hole masses for NGC 3842 and NGC 4889, 
or perhaps with their velocity dispersion measurements, 
or if instead these values are fine and these galaxies are 
representative (of a possibly new class) of galaxies which 
deviate from the M^-a relation as we know it, in the 
following section we have elected to repeat the M^-a 
regression analysis without these two galaxies. If these 
galaxies are outliers, then this additional approach pro- 
vides a more robust statistical analysis of the M^-a re- 



lation. To avoid issues of inconsistent samples, these two 
galaxies are also excluded from the ensuing M^-L dia- 
gram, although doing so has no significant effect on our 
M^-L relation. 

Although NGC 1316 is not an outlier from the M^-a 
diagram (Figure [lj , we do additionally exclude it from 
what follows due to the previously mentioned uncertainty 
associated with its core type. If NGC 1316 is confirmed 
to possess a partially-depleted core, something associated 
with dry merger events, then it will be the only such ex- 
ample in a barred galaxy (with bars of course associated 
with secular evolution rather than dry mergers). Until a 
three-component Sersic-bulge plus exponential-disk plus 
bar model has been fit to this galaxy, it is premature to 
claim from a fitted single-component, core-Sersic model 
that this galaxy has a depleted core. 

Our cautious exclusion of these three galaxies 
(NGC 1316, NGC 3842 and NGC 4889), leaves us with 
our final sample of (75-3=) 72 galaxies, comprised of 24 
core-Sersic galaxies and 48 Sersic galaxies. 

3.1. The Mbh-a diagram 

In Figure [2^, we have performed a symmetrical linear 
regression for the 24 core-Sersic galaxies in the M^-a 
diagram. This was achieved using the bisector regression 
analysis in the BCES routine from Akritas & Bershady 
(1996). The results are provided in Tabled] Of note 
is that the symmetrical M^-a relation has a slope con- 
sistent with a value of 5 (Ferrarese & Merritt 2000; Hu 
2008; Graham et al. 2011) . 

We have additionally performed the same symmetrical 
regression for the Sersic sample (see Table However, 
due to each supermassive black hole's limited sphere-of- 
influence, r- m a ~ GM hb /a 2 (Merritt & Ferrarese 2001b), 
there is an associated galaxy distance limit to which we 
can spatially resolve the black hole's dynamical influence 
on its surrounding gas and stars. Shown in Figured is 
the limit to which one can detect black holes at a dis- 
tance of 1.3 Mpc when the spatial resolution is 0".l (i.e. 
roughly that achievable with the Hubble Space Telescope 
and ground-based adaptive-optics-assisted observations) . 
For galaxies further away, this limit moves vertically in 
the diagram to higher black hole masses. Attempts to 
occupy the lower right of Figure [2^i may need to wait 
for the next generation of telescopes with greater spa- 
tial resolution, and as such we currently have a selection 
limit, a kind of floor, in this diagram (Batcheldor 2010; 
Schulze & Wisotzki 2011). Galaxies can exist below this 
floor, but they are at distances such that we can not re- 
solve their black hole's sphere-of-influence. To avoid this 
sample selection bias on the samples which extend into 
the lower portion of the diagram, we have additionally 
performed an ordinary least squares (OLS) regression of 
the abscissa (i.e., the horizontal X value) on the ordinate 
(i.e., the vertical Y value) for the reasons expounded in 
Lynden-Bell et al. (1988, their Figure 10). 

Given that the 27 non-barred Sersic galaxies define an 
Mtfo-a relation which is consistent with that defined by 
the 24 core-Sersic galaxies (see Table [3|), they have been 
combined to produce a single non-barred M^-cr relation. 
These 51 galaxies yield a slope of 5.53±0.34, consistent 
with the value of 5.32±0.49 reported by Graham et al. 
(2011) using a sample of 44 non-barred galaxies. Barth, 
Greene, & Ho (2005, their figure 2) offers further support 
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Fig. 1. — Black hole mass versus central velocity dispersion. The faint dot-dashed gray line corresponds to a black hole's sphere-of- 
influence of 0".l at a distance of 1.3 Mpc. With such a limiting spatial resolution of 0".l, black holes residing below this line cannot 
be reliably detected. Panel a) Red dots represent the core-Sersic galaxies, open blue circles represent the Sersic galaxies, while the open 
crosses designate those which are barred. All 77 galaxies shown here have been included in the analysis of this Figure. The upper solid red 
line corresponds to the symmetrical regression for the 27 core-Sersic galaxies, while (due to the sample selection limit) the dotted and solid 
blue lines correspond to the 'ordinary least squares' OLS(c|Mt,h) regression for the barred and non-barred Sersic galaxies (see Table [3j - 
Panel b) Red dots represent the non-barred galaxies while the open blue crosses designate those which are barred. The solid red line is a 
fit to the non-barred galaxies, while the dashed black line is a fit to all galaxies. M32, the Milky Way and Fornax A (NGC 1316) have been 
marked with a green square, they are not incl uded in Figure [2] for the reasons mentioned in Section [2] The Sersic galaxy with the highest 
black hole mass is NGC 3115 (see section . 



for this linear relation after reducing their active galactic 
nuclei (AGN) black hole mass estimates by a factor of 
^2 due to a revision in the virial /-factor used to deter- 
mine black hole masses in AGN (Onken et al. 2004; Gra- 
ham et al. 2011). While Graham et al. (2011) discussed 
why a 5% uncertainty on the velocity dispersion may be 
optimistic, we note that using a 10% uncertainty on er, 
rather than 5%, results in the same slope to the M^-a 
relation when using the OLS(cr|Mbh) regression. Cou- 
pling the knowledge that L oc a 2 for the Sersic spheroids, 
with the relation Mbh oc cr 5 5 , one would expect to find 
Mbh L 2 - 75 for the Sersic spheroids. 

For a value of a = 400 km s _1 , the non-barred M^-a 
relation yields a black hole mass of 7.7 x 10 9 M Q . This is 
3 times greater than the value of 2.6 x 10 9 Af Q predicted 
by the elliptical galaxy Mbh - cr relation from Giiltekin et 
al. (2009b), and better, although not fully, matches the 
expectations from Hlavacek-Larrondo et al. (2012) for 
black holes in brightest cluster galaxies if they reside on 
the "fundamental plane of black hole activity" (Merloni 
et al. 2003; Falcke et al. 2004). 

As reported by Graham (2007b, 2008a,b) and Hu 
(2008), the barred galaxies are, on average, offset to lower 
black hole masses or higher velocity dispersions than 
the non-barred galaxies, with the latter scenario possi- 
bly expected from barred galaxy dynamics (e.g. Gadotti 
& Kauffmann 2009; Graham et al. 2011; DeBuhr, Ma 
& White 2012). Consistent with Graham et al. (2011), 
we find a mean vertical offset of 0.30 dex between the 
barred and non-barred galaxies in the Mbh - c diagram, 
which corresponds to a factor of ~ 2 in black hole mass, 
when a — 200 km s _1 . 



For those who may not know if their galaxy of in- 
terest is barred or not, using the full 72 galaxies and 
an OLS(Mbh|c), bisector, and OLS(<r|Mbh) regression 
to construct the classical (all galaxy type) M^-a rela- 
tion gives slopes of 5.21 ± 0.27, 5.61 ± 0.27 and 6.08 ± 
0.31, respectively. The associated intercepts are 8.14 ± 
0.04, 8. 14 ± 0.05 and 8.15±0.05. These steeper relations 
are expected given the offset nature of the barred galax- 
ies. 

3.2. The M\fo-L diagram 

In Figures [2J5 and [2}: we have performed a symmet- 
rical linear regression for the 24 core-Sersic galaxies in 
the Muh-L diagram. Expressing their spheroid magni- 

11 n 1 n "7n 

tudes as luminosities, we find that iV/bh oc L ^ and 
Mbh oc £^ 35±0 - 30 ; j n reasonable agreement with the re- 
lation Mbh oc L 10 and with past analyses of predomi- 
nantly bright galaxies 10 . We note that depending on the 
progenitor galaxy mass ratios — among the individual 
dry merger events which built one's core-Sersic sample 
- the slope of the M^-L relation can be expected to 
deviate slightly from a value of 1. Although, for galaxies 
built from a sufficient number of dry mergers, this should 
become a second order effect in terms of the overall evo- 
lutionary scenario in the M^-L diagram, for the reason 
described by Peng (2007) and Jahnke & Maccio (2011). 

The effective sample selection boundary shown in Fig- 
ure [2^i has been mapped into Figures [2Jd and [2J: by con- 
verting this line's velocity dispersion into the expected 

10 The inclusion of NGC 3842 and NGC 4889 does not signifi- 
cantly change these results. 
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TABLE 3 
Black hole scaling relations 



# and Type Regression a /3 A dex 

Figure [l] (includes some questionable data points) 
log[M bh /M ] = a + /3 log[cr/200 kms" 1 ] 

27 Core-Sersic Bisector 8.02 ±0.17 6.95 ± 1.17 0.45 

28 Non-barred Sersic Bisector 8.21 ± 0.06 4.39 ± 0.53 0.32 
23 Barred Bisector 7.78 ±0.10 4.14 ±0.55 0.34 
54 Non-barred All Bisector 8.29 ± 0.05 5.09 ± 0.44 0.39 

28 Non-barred Sersic OLS(cr|M bh ) 8.25 ± 0.06 4.82 ± 0.64 0.34 

23 Barred OLS(cr|M bh ) 7.89 ±0.16 4.98 ± 1.00 0.37 
54 Non-barred All OLS(o-|M bh ) 8.28 ± 0.06 5.57 ±0.52 0.41 
77 All OLS(ojM bh ) 8.20 ±0.05 6.08 ± 0.41 0.46 

Figure [2] panel a), M bb -cr 
log[M bh /M Q ] = a + /3 log[ff/200 kms" 1 ] 

24 Core-Sersic Bisector 8.20 ±0.13 5.35 ± 0.92 0.35 

27 Non-barred Sersic Bisector 8.23 ± 0.06 5.02 ± 0.41 0.31 
21 Barred Sersic Bisector 7.75 ±0.12 4.03 ± 0.69 0.34 
51 Non-barred All Bisector 8.23 ± 0.05 5.14 ±0.31 0.33 

28 Elliptical Bisector 8.23 ± 0.05 4.86 ± 0.50 0.34 

27 Non-barred Sersic OLS(cr|Af bh ) 8.27 ±0.06 5.53 ± 0.51 0.33 
21 Barred Sersic OLS(<r|M bh ) 7.92±0.23 5.29±1.47 0.40 
51 Non-barred All OLS(<r|M bh ) 8.22±0.05 5.53±0.34 0.34 

28 Elliptical OLS(cr|M bh ) 8.21 ±0.10 5.54 ± 0.71 0.37 
72 All OLS(ojM bh ) 8.15 ±0.05 6.08 ± 0.31 0.41 

Figure[2] panel b), M bh -L K ^ aph 
log[M bh /M ] = <x + 0[M Ka:3ph + 25] 

24 Core-Sersic Bisector 9.05±0.09 -0.44±0.08 0.44 

log[M bh /M Q ] = a + P[M Ka>sph + 22.5] 

48 Sersic Bisector 7.39±0.14 -1.09±0.22 0.95 

Figurc[2] panel c), M bh -L BtSph 
log[M bh /M Q ] = a + i3[M Bt3ph + 21] 

24 Core-Sersic Bisector 9.03±0.09 -0.54±0.12 0.44 

log[M bh /M ] = a + P[M Bt3ph + 19] 

48 Sersic Bisector 7.37±0.15 -0.94±0.16 0.93 

M bb = black hole mass, a = host galaxy velocity dispersion, Mfc S} sph = dust-corrected, 2MASS A" s -band spheroid magnitude, Mg gp^ = 
dust-corrected, RC3 B-band spheroid magnitude. A symmetrical "bisector" regression was used, in addition to a non-symmetrical 
'ordinary least squares' OLS(A|Af bb ) regression to compensate for the sample selection limit, or floor in the data, at the low-mass end 

(see Figure [2}. The total rms scatter in the logM bb direction is given by A, and roughly scales linearly with the slope (within each 
panel). The lines highlighted in bold show the preferred fits, see the text for details, while Table [4] provides revised M bb -Lx s relations 
according to Schombert & Smith's (2012) correction of the 2MASS photometry. 



magnitude according to the L-a relations defined by 
the Sersic galaxies in these diagrams. The need for 
an OLS(L|Afbh) regression on the Sersic galaxies in the 
Mhh—L diagrams, rather than a symmetrical regression, 
is no longer as apparent because the 'floor' to the sample 
selection has significantly shifted/rotated. 

It is not clear if the barred Sersic galaxies are offset 
from the non-barred Sersic galaxies in Figures [2b and 
[2J; and so we have therefore grouped them together in 
our analysis. The results of which, given in Table [3j are 
such that M bh oc L 2 K ] 3±0 - 55 and M bh oc L 2 B 35±0A0 for the 
Sersic spheroids 11 . We note in passing that it is expected 
that the color-magnitude relation for Sersic spheroids 
(e.g., Tremonti et al. 2004; Jimenez et al. 2011), as op- 
posed to the flat color-magnitude relation for core-Sersic 
galaxies (as discussed in Graham 2012c), will result in 
these two exponents not being equal to each other. How- 
ever the current uncertainty on these two exponents does 
not allow us to detect this difference. 

Finally, Schombert (2011) and Schombert & Smith 
(2012, see their Appendix) have recently detailed a pho- 

11 The inclusion of NGC 3079 does not change these results. 



tometry error in the 2MASS Extended Source Catalog 
(Jarrett et al. 2000), such that the total 2MASS magni- 
tudes are, on average, 0.33 mag too faint in the J-band. 
There is a similar offset in the -£T s -band data (Schombert 
2012, priv. comm.) such that the V — K s colors are too 
blue by this amount. Therefore, for those using K or 
-FTg-band data not from the 2MASS catalog for the pre- 
diction of black hole masses, the normalizing values of 
25 and 22.5 used in the M^-Lks equations shown in 
Table [3] should be increased by 0.33 to 25.33 and 22.83. 
One can actually go a little further than this. Assigning a 
1-sigma uncertainty of 0.33/2 — 0.165 mag to this correc- 
tion, that is, assuming that the scatter seen in Schombert 
& Smith's figure 12 has a 1-sigma value of 0.33/2 mag, 
we can estimate the impact of this correction on our K s - 
band M^-L relations. Because they were constructed 
using N — 24 core-Sersic galaxies and N = 48 Sersic 
galaxies, our 'normalizing' values of 25 and 22.5 mag have 
an additional uncertainty of 0.165/a/24 and 0.165/-\/48 
mag associated with them, which is a small but non- 
zero 0.034 and 0.024 mag, respectively. Given the slopes 
of our two Mbh _ iif 3 relations, this is equivalent to in- 
creasing the uncertainty on their intercept from 0.09 to 
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TABLE 4 

Corrected M hh -Lx s scaling relations 



# and Type Regression 


a 





A dcx 


log[A/ bh /A%] = a + 
24 Core-Sersic Bisector 


P[Mli s ,sph ■ 

9.05 ± 0.12 


+- 25.33] 
-0.44 ±0.08 


0.44 


log[M bh /A%] = a + 
48 Sersic Bisector 


P[^K a ,sph ■ 

7.39 ± 0.16 


■f 22.83] 
-1.09 ± 0.22 


0.95 



X s -band M^-L relations after correcting for the 2MASS 
magnitude errors noted by Schombert & Smith (2012). 

0. 12 and from 0.14 to 0.16, respectively. The result of 
applying these corrections is shown in Table 0| . While 
this is a rather minor adjustment to our equations, it 
remains true that those using the 2MASS catalog of K s - 
band magnitudes for individual galaxies may have the 
true galaxy magnitude wrong by the extent shown by 
Schombert & Smith (2012), which will thus affect their 
predicted black hole mass. 

It may be worth providing one additional comment 
in regard to photometric errors. If the 2MASS mag- 
nitudes are shown to contain a systematic error such 
that a greater fraction of galaxy light at large radii is 
increasingly missed in the intrinsically brighter galaxies, 

1. e. those with extended light profiles, this will not ex- 
plain the bend in the M^-L relation. Bulges in (often 
truncated) exponential disk are not significantly affected 
by this potential problem. The affect is greater in galax- 
ies with higher Sersic indices, specifically the core- Sersic 
elliptical galaxies. Such a correction would therefore act 
to slightly reduce the slope of the core-Sersic M^-L re- 
lation, thus increasing the apparent 'bend' in the M^—L 
diagram between the Sersic and the core-Sersic spheroids. 

3.3. Propagation of errors 

It is important to appreciate the uncertainties on the 
slope and intercept of the relations given in Tables [3] and 
IH plus the associated intrinsic scatter e. For example, the 
maximum 1-sigma uncertainty on the predicted value of 
Mbh using the M^-Lk, relation is acquired by assuming 
uncorrelated errors on the magnitude Mk s and the slope 
and intercept of the relation. Gaussian error propagation 
for the linear equation y = (b ± 5b) (x ± 5x) + (a ± <5a), 
gives an error on y equal to 

5y= y / (dy/db) 2 (5b) 2 + (dy/da) 2 (5a) 2 + (dy/dx) 2 (Sx) 2 

= v / x 2 {5b) 2 + (5a) 2 +b 2 (5x) 2 . 

In the presence of intrinsic scatter (e) in the y-direction, 
the uncertainty on y is 

8y = v /x 2 (Sb) 2 + (5a) 2 + b 2 (5x) 2 + e 2 . 

For the Sersic Mhh~LK s relation in Table SJ x = 
[M Ka>ap h + 22.83], and thus' 

(5 log M hh /M Q ) 2 = [M KatSph + 22.83] 2 (0.22) 2 + (0.16) 2 

+(-1.09) 2 [SM Ks , sp h} 2 + e 2 , 

where 8MK 3 , S ph is the uncertainty associated with the 
spheroid's magnitude. 

For the preferred non-barred M^-a relation in Ta- 
ble H x = log(cr/200kms" 1 ), so dx/da = l/[ln(10)cr], 
and therefore 

(<51ogAf bh /M ) 2 = [log(cr/2OOkms" 1 )] 2 (O.34) 2 + (0.05) 2 
+ [5.53/ln(10)] 2 [5cT/a] 2 + (e) 2 . 



Precise knowledge of the intrinsic scatter in the above 
equations is a difficult task because it requires an accu- 
rate knowledge of the measurement errors, which we do 
not claim to have. As such, the scatter in our current 
Mbh-L diagram probably cannot be used to help dis- 
tinguish between competing formation scenarios, as has 
been cleverly proposed (e.g. Lahav et al. 2011; Shankar et 
al. 2012). The measurement errors should add in quadra- 
ture with the intrinsic scatter to give the total rms scatter 
A. We can thus quickly proceed with a rough estimate 
of the intrinsic scatter. In terms of the vertical scatter 
about the M^-a relation, we have seen from Graham et 
al. (2011, their Table 2) that e « 3A/4. Here we shall 
use a rough value of 0.3 dex, which is consistent with 
the results in Gultekin et al. (2009b). For the M hh -L 
relation, Graham (2007a, his Table 4) reported a similar 
value of around 0.3 dex for his predominantly luminous 
galaxy sample. Obviously studies which have attempted 
to measure this (see also Gultekin et al. 2009b) but failed 
to account for the bent nature of the M^—L relation will 
be in error at some level. We shall however use this value 
here for the core-Sersic galaxies. For the fainter Sersic 
galaxies, which follow an M^-L relation that may be 
some 2.5 times steeper, we use a value of 0.75 dex. This 
increase to the intrinsic scatter in the vertical direction 
is expected if the intrinsic scatter in the horizontal direc- 
tion, i.e. for the magnitudes, is equal for both types of 
galaxy, i.e. core-Sersic and Sersic galaxies. 

The above equations are used in Section [4 . 1 1 where we 
predict a number of black hole masses which are less than 
one million solar masses. 

3.4. Method of regression 

Park et al. (2012) explored how the BCES linear re- 
gression code from Akritas & Bershady (1996) can pro- 
duce a biased M^-a relation if one's velocity dispersions 
have large measurement errors. They revealed that there 
is not a bias when these errors are less than ~15%, the 
uppermost value used in the literature. Park et al. (2012) 
additionally showed that the 'forward' and 'inverse' lin- 
ear regression from the BCES code yield M^-a rela- 
tions which are consistent with those derived using the 
'forward' and 'inverse' versions of the modified FITEXY 
expression from Tremaine et al. (2002). That there is 
agreement between these codes in the M^-a diagram is 
not a new result (e.g. Novak et al. 2006; Graham & Li 
2009; Graham et al. 2011), but it may be a poorly appre- 
ciated point worthy of some explanation before using the 
modified FITEXY code to check on our BCES-derived 
M^h~L relations. 

It is known that a linear regression which minimizes 
the offset of data from a fitted line will typically yield 
a different result if one minimizes the offset in either 
the vertical or the horizontal direction (e.g. Feigelson & 
Babu 1992). It is also known that minimization of the 
vertical offsets (the so-called 'forward' regression) yields 
a shallower slope than obtained through minimization 
of the horizontal offsets (the 'inverse' regression). Ob- 
viously neither of these are symmetrical regressions be- 
cause there is a preferred variable, that is, they do not 
treat the x and y data equally and therefore they gen- 
erate different fitted lines. The two lines from each of 
these non-symmetrical regressions can however be com- 
bined to produce a symmetrical regression. For example, 
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Fig. 2. — Black hole scaling relations involving a) central velocity dispersion a, b) i\s-band host spheroid magnitude and c) S-band host 
spheroid magnitude. The red dots represent the (luminous) core-Sersic galaxies and the blue circles represent the (intermediate-luminosity) 
Sersic galaxies and bulges, while the open crosses designate those which are barred. The symmetrical bisector regressions given in Table [3] 
are shown here, except for panel a) which shows the OLS(<r| Afth) regression for the 21 barred (dotted blue) and 27 non-barred (solid blue) 
Sersic galaxies, along with the symmetrical bisector regression for the 24 core-Sersic galaxies. The faint dot-dashed gray line at the bottom 
of panel a) corresponds to a black hole's sphere-of-influence of 0".l at a distance of 1.3 Mpc. With such a limiting spatial resolution of 0".l, 
black holes residing below this line cannot be reliably detected. Using the Lk s -& and L^-cr relations for the Sersic galaxies, constructed 
using both a bisector regression and an OLS(L|ct) regression, this gray line has been mapped into panels b) and c). The left-most core-Sersic 
galaxy in panels b) and c) is NGC 3998. 



the line which bisects the forward and inverse regression 
lines, having the average angle between them, is the sym- 
metrical bisector regression line. In general, symmetri- 
cal regressions are preferred when it is not known which 
variable is dependent on the other, and when one wants 
to establish the intrinsic relation between two quanti- 
ties, the so-called theorist's question (Novak et al. 2006). 
Biases in the data sample, due to selection boundaries, 
can of course bias the results from a symmetrical regres- 
sion and thus create a preference for either the forward 
or inverse regression depending on whether the selection 
boundary is vertical or horizontal in one's diagram (e.g. 
Lynden Bell et al. 1988). Typically in the literature, the 
modified FITEXY routine has only been run as a 'for- 
ward' regression, minimizing the offsets in the vertical 
direction, where as users of the BCES code typically re- 
port the symmetrical bisector regression line. It is not 
appropriate to compare such results, and some confu- 
sion exists because Tremaine et al. (2002) remarked that 
their modified FITEXY routine treats the data symmet- 
rically. However if this was true, if the x and y data was 
treated equally, then swapping the x and y data points 
would not alter the fit that one obtains. 

Returning to our M^-L diagram, given the larger un- 
certainties on the spheroid magnitudes, compared to the 
uncertainties on the velocity dispersions, we have used 
the modified FITEXY routine to check if the BCES 
bisector routine has yielded a biased result. In fitting 
the line y = a + bx, Tremaine et al.'s (2002) modified 
version of the routine FITEXY (Press et al. 1992, their 
Section 15.3) minimizes the quantity 
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bxi 
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x? + e 2 



(8) 



The intrinsic scatter (in the y-direction) is denoted here 
by the term e, and the measurement errors on the ./V 



pairs of observables Xi and are denoted by Sxi and 
Syi. The intrinsic scatter e is solved for by repeating the 
fit with different values of e until X 2 /(-/V — 2) equals 1. 
To achieve a minimization in the rr-direction, i.e. to per- 
form the 'inverse' regression, one can simply replace the 
e 2 term in the denominator with b 2 e 2 (Novak et al. 2006; 
Graham 2007a; Graham & Driver 2007a). The M bh -L 
bisector line for the core-Sersic spheroids, obtained from 
the bisector line of the forward and inverse modified FI- 
TEXY routine, has a slope of —0.46 and —0.56 in the 
K- and B- bands, respectively. The slope of the M^-L 
bisector line for the Sersic spheroids, obtained from the 
forward and inverse modified FITEXY routine, is —0.95 
and —0.88 for the K and -B-bands, respectively. These 
results are consistent with the values obtained from the 
BCES bisector routine, as reported in Table [3] 

3.5. Further Pruning 

Curious readers may be wondering what happens if 
additional data points are pruned from the M^-L dia- 
gram. The 'core' galaxy with the faintest magnitude is 
the lenticular galaxy NGC 3998. With a velocity disper- 
sion in excess of 300 km s _1 , this galaxy is expected to 
be a 'core' galaxy. Figure 8 from Gonzalez Delgado et 
al. (2008) reveals that this galaxy hosts an AGN that 
swamps the central light profile and any core which may 
be there, hence the cautious '?' on its classification in 
Table 1. Re-labelling this galaxy to be a "Sersic" galaxy 
slightly steepens the Sersic M^-L relation, and hence 
our case for a deviation from the previous linear M-L re- 
lation, but not significantly. To the immediate lower left 
of this galaxy in the M hh -L diagram is NGC 4342. This 
Sersic galaxy was labelled an outlier by Laor (2001) and a 
potential outlier by Haring & Rix (2004), both of whom 
excluded it from their analysis. Bogdan et al. (2012a) 
have also identified it as unusual, and it may be a heavily 
stripped galaxy (Blom et al. 2012), causing it to deviate 
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Fig. 3. — Schematic view of Figure[2p showing the evolutionary 
path of 'dry' galaxy mergers as they branch off from the steeper, 
near-quadratic, M^—L relation for Sersic galaxies built from 'wet' 
mergers and/or gas-rich processes. The upper red arrow roughly 
shows the direction of growth if oc L 10 , while the slope of the 
blue shaded area roughly follows Mbh x L 2 ' 5 - An additional five 
galaxies (NGC 3842 and NGC 4889 from McConnell et al. 2011a, 
NGC 1316, NGC 3079 and M32), which were excluded from the 
regression analysis in Figure [2] are shown here, along with O J278 
from Valtonen et al. (2012), and marked with a green square. 

from the M^-L relation while still following the Mbh~ 
a relation. M32 may be another such stripped galaxy 
(Bckki 2001; Graham 2002). Excluding NGC 3998 and 
NGC 4342 results in an M\^-Lk b relation with slopes of 
-0.50 ± 0.08 and -1.07 ± 0.24 for the Sersic and core- 
Sersic galaxies, respectively. 

4. DISCUSSION AND IMPLICATIONS 

The coevolution of black holes and their host spheroids 
has resulted in a bent, i.e. a non-(log-linear), M^-L re- 
lation which can be approximated by a broken power-law 
having an exponent of ~1 for spheroids with core-Sersic 
profiles and ~2.5 for spheroids with Sersic profiles. The 
major, i.e. near equal mass, merger of gas-free galax- 
ies near the bright-end of the Sersic sequence will pre- 
serve the Mbh/L ratio and create the core-Sersic galax- 
ies, which are observed to follow a near linear M^-L 
relation. That is, such dry mergers do not move galaxies 
along the steep M^-L relation defined by Sersic galax- 
ies, but drive them off it to form a distinctly different 
sequence (Figured]) 12 . The above picture (and Figured]) 
suggests that the bulk of AGN emission, i.e. black hole 
growth through gas-rich processes, should be associated 
with black holes having masses less than a few times 
1O 8 M . The peak of the Eddington ratio at M hh =(4- 
8)xl0 7 M Q is not inconsistent with this (DeGraf et al. 
2012). Magnitude-limited high-redshift quasar surveys 
may thus only be probing the tip of the proverbial ice- 

12 A substantial gas-rich merger with, or gas inflow onto, a core- 
Sersic galaxy may drive it above the one-to-one relation seen in the 
M^h-L diagram. 



berg. 

An exponent of 1 for the core-Sersic galaxies in the 
M\^-L diagram would mean that there is a constant 
Afbh/i^spheroid mass ratio for these galaxies. For a core- 
Sersic spheroid with a (dust-free) absolute magnitude 
Mk, = —25.33 mag, and assuming an M/Lk ratio of 0.8 
(Bell & de Jong 2001) and a solar if -band absolute mag- 
nitude of 3.27 mag (Cox 2000), we find a median mass 
ratio Mbh/Msph,* = 0.49%. For a core-Sersic spheroid 
with a (dust-free) absolute magnitude Mb = — 21 mag, 
and assuming an M/Lb ratio of 5 and a solar B-band 
absolute magnitude of 5.47 mag (Cox 2000), we find 
a similar (black hole)-to- (stellar) mass ratio of 0.55%. 
These ratios are larger than the old median value of 0.14- 
0.2% (e.g. Ho 1999; Marconi & Hunt 2003; Haring k 
Rix 2004) for the (black hole)-to- (spheroid dynamical) 
mass ratio. Even if we assumed M/Lb = 8, we find 
M bh /M Sph ,* = 0.35%. 

The steep exponent of ^2.5 for the Sersic galaxies in 
the M\,\y-L diagram means that the pre-existing single 
log-linear M^—L relation will overestimate the black hole 
masses in galaxies whose black hole mass is less than 
^10 8 Af Q . This is in fact already evident in Khorun- 
zhev et al. (2012) who found that Mbh was overestimated 
using the old M^-L relations for Mbh < 10 8 Af Q , just 
where the bend in the M^—L relation sets in. Similarly, 
Mathur et al. (2012, their figure 3) show ten galaxies be- 
low the old Mbh-i relation at log(Mbh) = 7± 0.4. How- 
ever their interpretation that these galaxies host pseu- 
dobulges, based on their location in the M^-L diagram, 
is premature because these galaxies are found where the 
classical bulges and elliptical galaxies reside. 

Having the correct black hole scaling relations is, obvi- 
ously, important when trying to predict a galaxy's black 
hole mass. Equally important is having the correct ve- 
locity dispersion or bulge magnitude. In regard to the 
core-Sersic galaxy NGC 4382 (Dullo & Graham 2012), 
Giiltckin et al. (2011) predicted a black hole mass that 
differed by an order of magnitude when one used the 
M^-L and Mbh-c relations (Giiltekin et al. 2009b). 
Here we predict log(Mbh/M©) = 8.0 ± 0.3 when using 
a = 182 ± 5 km s _1 (Giiltekin et al. 2011), in agree- 
ment with Giiltekin et al. (2011). However this galaxy 
appears to be a disturbed, face-on, non-barred, lentic- 
ular galaxy (Laurikainen et al. 2011) and we therefore 
do not use the total galaxy light to estimate its central 
black hole mass. The total 2MASS if s -band magnitude 
of 6.145 mag corresponds to an absolute magnitude of 
roughly —25.1 mag given a distance modulus of 31.27 
(Tonry et al. 2001; Blakeslee et al. 2002). Galactic ex- 
tinction and redshift dimming are negligible. Using an 
average, dust-corrected, SO galaxy bulge-to-disk flux ra- 
tio of one third, we thus have a if s -band bulge magnitude 
of —23.9 mag. Using the M^-Lks relation in Table [3j 
we find log(Aibh/-M©) = 8.6 ± 0.5 when using a factor of 
2 uncertainty on this bulge magnitude. This estimated 
black hole mass has over-lapping error bars with the first 
estimate. We would however predict a black hole mass 
of ~ 10 9 A/ Q , in agreement with Giiltekin et al. (2011), if 
all of the light in this disturbed galaxy was assigned to 
a single spheroid component. 

Having the right M^-L relation is of course also im- 
portant when using the galaxy/spheroid luminosity func- 
tion to predict the black hole mass function, and sub- 
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sequently the black hole mass density (e.g. Graham & 
Driver 2007b, and references therein). Knowledge of 
the correct M^-L relation will additionally be useful 
for cosmological / evolutionary studies of the nucleus-to- 
(host spheroid) mass ratio (e.g. Kisaka & Kojima 2010; 
Schulze & Wisotzki 2011; Portinari et al. 2012; Sesana 
2012), and for many other studies. The consequences of 
the new bent M^-L and bent Mbb-Mgphcroid relations 
are widespread and dramatic. Below we discuss some of 
these in more detail. 

4.1. Intermediate mass black holes 

Classical elliptical galaxies with Mb ~ —18 mag will 
now have (Mbh _ £)-derived black hole masses which are 
ten times lower than previously predicted. Similarly, 
we also now have revised predictions for the location of 
10 6 M© black holes within the M^-L diagram. Based 
on their host spheroid luminosity, they will have black 
hole masses which are roughly 10 times lower than pre- 
dicted from the extrapolation of the M^-L relation de- 
fined by luminous core-Sersic galaxies, placing them near 
the current detection limit shown in Figurc[21 Black holes 
with masses of 10 5 -10 6 Mq will be even further from the 
extrapolated core-Sersic M^-L relation. A number of 
AGN residing in elliptical galaxies, not to be confused 
with the bulges of disk galaxies, have been observed be- 
low this relation. For example, Greene et al. (2008) 
showed some AGN which had black hole masses that 
are an order of magnitude below both the old M^-L 
relation and the Mbh _ Af S phcroid relation from Haring & 
Rix (2004). Our interpretation of this offset differs from 
Greene et al. (2008), in that we consider their sample 
of non-disk galaxies to be the extension of the Sersic se- 
quence, rather than a separate galaxy population which 
Greene et al. (2008) refer to as "blobs" . 

4.1.1. Individual galaxies 

Interestingly, the (Mbh-Af S pheroid)-derived mass deter- 
mination (using the expression from Haring & Rix 2004) 
for the black hole associated with the recent ultraviolet- 
optical flare from the tidal disruption of a helium-rich 
stellar core in a galaxy at a redshift of 0.17 (Gezari 
et al. 2012) can now be revised. Using the updated 
-Wbh-^sphcroid relation from Graham (2012a), its black 
hole mass is ~5 times lower than the reported value of 
4 x 10 6 M Q , making it slightly less than one million solar 
masses. 

Our revised black hole scaling relations suggest that 
the high-energy transient Sw 1644+57 (Levan et al. 2011; 
Bloom et al. 2011) may have been associated with a 
stellar capture event by a black hole having a mass of 
0.6 x 10 6 M Q (and even less if the host galaxy is a disk 
galaxy). Not only does the host galaxy's absolute mag- 
nitude (My = —18.19 mag) suggest this low black hole 
mass (for B — V — 0.9), but so does the reported host 
galaxy mass of a few billion solar masses (Bloom et al. 
2011) when using the -Mbh--M sp heroid relation from Gra- 
ham (2012a). We therefore tentatively claim the presence 
of an intermediate mass black hole in this galaxy, adding 
to the detection of other intermediate mass black hole 
candidates (e.g., Farrell et al. 2009; Soria et al. 2010; 
Sutton et al. 2012). 

The revised M^-L scaling relations presented here 
give rise to the need to re-derive many past estimates of 



black hole masses, in particular those associated with ac- 
tive galactic nuclei in low luminosity bulges. Whether or 
not these objects can provide unique observational con- 
straints on the progenitors of supermassive black holes is 
unclear, but they do narrow the divide between today's 
supermassive black holes and their 'seeds'. 

For the Seyfert Sc galaxy NGC 3367 (McAlpine et al. 
2011), Dong & De Robertis (2006) report a if -band bulge 
magnitude of —21.6 mag, corresponding to a bulge-to- 
total flux ratio of 0.08. Hernandez-Toledo et al. (2011) 
find that the bulge has a Sersic index between 1 and 
1.8, and after including a bar in their image analysis 
they report a bulge-to-disc ratio of 0.07 to 0.11. The 
bulge magnitude from Dong & De Robertis (2006), after 
brightening it by 0.11 mag for dust (Driver et al. 2008), 
yields an expected black hole mass of log M^/Mq — 
6.2 ± 0.9. McAlpine et al. (2011) also studied the Sbc 
galaxy NGC 4536, reported to have a bulge Sersic index 
of 1.88 ±0.35 by Fisher & Drory (2008), and an absolute 
V-band magnitude of -18.4 mag assuming a distance of 
25.3 Mpc. With a disk inclination of 62 degrees (Fisher 
& Drory 2008) we apply a dust correction of 1.3 mag 
(Driver et al. 2008), and using an assumed B — V bulge 
color of 0.8 we have that Mb = —18.9 mag. The ex- 
pected black hole mass is thus log-Mbh/Ms = 7.3 ± 0.8, 
assuming a 0.3 mag uncertainty on the bulge magnitude. 
These two galaxies are not, therefore, expected to house 
low intermediate mass black holes. 

If the dwarf galaxy Pox 52 — long thought to house an 
AGN (Kunth, Sargent & Bothun 1987) — is an elliptical 
galaxy, then its -B-band magnitude of —16.8 mag (Barth 
et al. 2004) translates into an expected black hole mass 
of logMbh/M© = 5.3 ± 0.9, with the uncertainty of 0.5 
dex derived assuming an uncertainty of 0.3 mag on the 
galaxy magnitude. 

For the nearby dwarf lenticular galaxy NGC 404, 
the Galactic extinction corrected 7J-band magnitude of 

— 19.65 mag from Scth et al. (2010) becomes Mk s — 

— 19.85 mag when using an H — K s color of 0.2 mag (e.g. 
Balcells et al. 2007). Thought to house an AGN (e.g. 
Nyland et al. 2012, and references therein), the expected 
black hole mass in NGC 404 is logM b h/M Q = 4.1 ± 1.1, 
using a magnitude error of 0.3 mag. For comparison, al- 
though restricted by assumptions about the inner stellar 
mass-to-light ratio, the stellar dynamical modeling from 
Seth et al. (2010) suggests that NGC 404 has a black 
hole mass less than 1O 5 M0. 

The bulge of the Scd galaxy NGC 1042 has a dust 
corrected magnitude of — 19.13 ± 0.75 K-m&g (Knapen 
et al. 2003; Graham & Worley 2008) and houses an active 
black hole (Shields et al. 2008; Seth et al. 2008). From 
the new if-band A/BH-^buigo relation, the expected black 
hole mass is log M^/Mq = 3.3 ± 1.4. 

Given a _B-band bulge magnitude of —13.93 mag 
(Graham & Spitler 2009) for the AGN-hosting galaxy 
NGC 4395 (Filippenko & Ho 2003), the expected black 
hole mass in this late- type galaxy's bulge is equally small. 
Assuming a 0.75 mag uncertainty on the bulge magni- 
tude, the black hole mass is log M^/Mq = 2.6± 1.3. Al- 
lowing for the 1-sigma uncertainty given here, this mass 
resides just below the 10 4 -10 5 M o mass range estimated 
by Filippenko & Ho (2003) and is also below the (rever- 
beration mapping )-derived estimate of 3.6 ± 2.2 x 1O 5 M 
(Peterson et al. 2005), even after it has been halved to 
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TABLE 5 

SUPERMASSIVE AND INTERMEDIATE MASS BLACK HOLES 



Ident. magnitude log Mbh 



NGC 4536 


-18.9B-mag 


7.3 


± 


0.8 


NGC 3367 


-21.71 K- 


-mag 


6.2 


± 


0.9 


Gezari et al. 






5.9 


± 


0.9 


Sw 1644+57* 


-18.19 V- 


mag 


5.8 


± 


0.9 


Pox 52 


-16.8B-mag 


5.3 


± 


0.9 


NGC 404 


-19.85 K- 


-mag 


1.1 


± 


1.1 


NGC 1042 


— 19.13 K- 


-mag 


3.3 


± 


1.4 


NGC 4395 


-13.93 B- 


mag 


2.6 


± 


1.3 



* This transient's designation is used to identify the host galaxy. 

1.8 ± 0.6 x 1O 5 M due to the updated virial /-factor 
in Graham et al. (2011). A lower black hole mass is 
however expected from the short time-scale variations of 
the X-ray flux which can vary by an order of magnitude 
(Vaughan et al. 2005), and would of course necessitate 
a higher Eddington accretion rate than the previously 
assumed maximum of just a couple of percent. 

All of our black hole mass estimates above have been 
collated in Table [5] for ease of reference and overview. 

Typically, when X-ray flux is detected in galaxies hav- 
ing a candidate intermediate mass black hole, the level of 
emission is so low that skeptics point to possible stellar 
origins for its production. That is, in general, there is 
not the danger that such galaxies might have too much 
X-ray luminosity supportive of a supermassive black hole 
rather than an intermediate mass black hole. 

4.1.2. A catalog of galaxies: Dong & De Robertis 

In addition to the above mentioned galaxies which have 
had individual papers dedicated to them, we are able to 
identify many other intermediate mass black hole can- 
didates from a catalog of low-luminosity AGN. The log- 
linear i^-band M^-L relation of Dong & De Robertis 
(2006), which was consistent with that from Marconi & 
Hunt (2003), was used by Dong & De Robertis to predict 
the masses of black holes in 117 disc galaxies hosting low- 
luminosity AGN. From their sample of bulge luminosi- 
ties, they reported only 4 galaxies with estimated black 
hole masses less than 10 6 M© and such that logMbh/M 
= 5.6 to 5.8. Given this, coupled with the uncertainty 
in the M^-L relation, they understandably paid no at- 
tention to the possible detection of intermediate mass 
black holes in their sample. However given the dramat- 
ically steeper M^—L relation reported in this paper for 
intermediate- and low-luminosity spheroids, we expect 
that many of the black hole mass estimates in Dong & 
De Robertis will have been overestimated. As a con- 
sequence, we note in passing that the associated scal- 
ing relations reported there between black hole mass and 
line width (FWHM [N II]), as well as emission-line ratios 
(e.g., [O III]/H/3, [O I]/Ha, [N II]/Ha, and [S II]/Ha), 
may need to be re-derived. 

Here we predict new black hole masses for the 41 galax- 
ies reported by Dong & De Robertis (2006, their Ta- 
ble 4) to have absolute if s -band bulge magnitudes fainter 
than —21.0 mag. Before applying our Sersic M^-L re- 
lation, we correct (brighten) these bulge magnitudes for 
internal dust extinction. This has been done using the 
inclination-dependent, dust correction from Driver et al. 
(2008), and is shown in Tabled The predicted black 
hole masses shown there assumes that the M^-L rela- 



tion can be extrapolated to lower black hole masses, and 
the quoted error was derived assuming an uncertainty on 
the magnitudes of 0.3 mag. 

The two lowest mass, intermediate mass black 
hole candidates in Table [B] are the barred galaxies 
NGC 4136 (log M bh /M Q = 2.2 ± 1.3) and NGC 3756 
(logMbh/M = 2.6 ± 1.3). We can apply our updated 
Mbh-c relation for barred galaxies, shown in bold in Ta- 
ble [21 although a serious concern is that the available 
velocity dispersions of 38.4±8.7 and 47.6±8.5 (Ho et al. 
2009) may be elevated by the velocity dispersion of the 
disk given the low bulge luminosities. The respective 
(Mbh-c)-derived masses logMbh/M — 4.13 ± 1.26 and 
4.62 ± 1.11 do however still have over-lapping error bars 
with the (Afbh~£)-derived masses. 

Finally, we note that moving down the spheroid lumi- 
nosity function, black holes appear to give way to, or 
at least co-exist with, dense nuclear star clusters (e.g. 
Graham & Spitler 2008; Graham 2012b; Neumayer & 
Walcher 2012; Leigh, Boker & Knigge 2012; Scott & 
Graham 2012). It is not clear to what lower masses our 
Sersic M^-L relation may hold. We do however note 
that many new intermediate mass black hole candidates, 
with Mbh = 0.8 x 10 5 -10 6 Af Q , measured from virial mass 
estimator based on the broad Ha line, have recently been 
identified (Dong et al. 2012). 

4.2. The Mbh-Afhost-gaiaxy connection 

For early-type galaxies, the dynamical mass-to-light 
ratio is thought to roughly change with the optical 
spheroid luminosity as (M/L)d yn amicai oc L 1 / 4 (e.g. Faber 
et al. 1987; Cappellari et al. 2006, their eq. 17) or rather 
cx L 1 / 3 (Cappellari et al. 2006, their eq.9 and 11), while 
MjL should obviously be constant for the core-Sersic 
galaxies if built via dry (additive) mergers. That is, the 
core-Sersic Mbh-Ms p h,dyn relation should have the same 
(close to unity) slope as the core-Sersic Mbh-£ sp hcroid re- 
lations Mbh oc L^ 35±0 - 30 and M hh oc £^ 0±a20 . However 
we need to caution that galaxy samples with a greater 
fraction of core-Sersic galaxies relative to Sersic galax- 
ies may have in the past resulted in a shallower average 
(M/L)cj y namicai oc L relation than is actually applicable 
to the Sersic galaxies on their own. Studies of the Fun- 
damental Plane (Djorgovski & Davis 1987) in the near- 
infrared J- and X-bands, which are less affected by stel- 
lar population differences across the galaxy luminosity 

range, have tended to yield (M/L) dynamical oc l][ 6 (e.g. 
Magoulas et al. 2012; La Barbera et al. 2010), although 
the expanded SAURON survey with 72 E/SO/Sa galaxies 

reported (M/L) dynamica i oc L\/ 3 and oc L\[^ (Falcon- 

Barroso et al. 2011). The ATLAS 3D survey (Cappel- 
lari et al. 2011) is comprised of just a few percent of 
core-Sersic galaxies supported by random motions rather 
than rotational velocity (Emsellem et al. 2011) and it 
will therefore be interesting to see if they also find these 
steeper M/L-L relations. 

For now, using the shallower MjL relations above, 
which admittedly may underestimate the slope of the 
Sersic galaxy (M/L)d yn amicai oc L relation, our Sersic 
Mbh _ ^spheroid expressions can be roughly transposed into 
(black hole)-(host spheroid) mass relations. One has 
that M bh oc L| 3 5±0-40 maps mt0 Mhh K M^ 8 ™ 2 
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TABLE 6 

Predicted intermediate mass black holes 
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Mx s ,orig is each galaxy's A" s -band, spheroid magnitude reported 
by Dong & De Robcrtis (2006). Mjf i)Corr is the dust-corrected 
magnitude following Driver et al. (2008). M^h is our predicted 
black hole mass using the Sersic M^-Lk 3 relation from Table [4] 

(For reference, Dong & De Robertis (2006) predicted only 4 
galaxies to have black hole masses less than 1O 6 M0, which were 
in the range log Af^h/ 'Mq = 5.6 to 5.8.) 

(when using M/L B oc L^ 4 ) and M bh oc L^3±o.55 

maps into Mbh oc M 2 ^ 4 ^ 47 (when using M/Lk b oc 

L^ 6 ). This can be directly compared with the slopes of 
1.01±0.52 for the core-Sersic galaxies in Graham (2012a), 
and 2.30±0.47 and 1.92±0.38 for the Sersic and non- 
barred Sersic galaxies in that study. 

Rather than described by a single power-law (e.g. 
Magorrian et al. 1998), the full Mbh-Ms p h,dyn relation 
appears to be better described by a broken power-law 
having exponents of ~1 for the core-Sersic spheroids and 
^2 for the Sersic spheroids, in agreement with the pre- 
vious announcement by Graham (2012a) based on a re- 
analysis of the Haring & Rix (2004) data (see also Bal- 
cells et al. 2007, their equation 13). The non-linear, 
albeit still a single log- linear, (black hole)- (host stellar 



mass) relation from Laor (2001) was given as Mbh oc 
^a.53±o.i4 an( ^ ^ mistake but rather arose from 

bph,* 

fitting a sample of what are undoubtedly Sersic and core- 
Sersic spheroids. 

If, as commonly believed, the core-Sersic galaxies are 
built from the simple addition of high-mass gas-free 
Sersic galaxies and/or other core-Sersic galaxies, then 
one may similarly expect Mbh to be linearly correlated 
with the dark halo mass in core-Sersic galaxies (Fer- 
rarese 2002; Bogdan et al. 2012a). Of course this (non- 
causal) relation will not hold across the full mass range 
of Sersic spheroids (Ho 2007). The single power-law re- 
lation involving the mass of the central black hole and 
the surrounding dark matter halo which was derived for 
a sample of both Sersic and core-Sersic galaxy types, 
Mbh oc M^ 5-1 - 82 (Ferrarese 2002), may be better de- 
scribed by a broken power-law — at least over the high- 
mass range where a correlation exists. Looking at Fig- 
ure 5 from Ferrarese (2002), galaxies with Mbh ^ 10 8 Af Q 
appear to be better described by a relation with an expo- 
nent closer to 1 than 1.65-1.82, with a quick analysis of 
their Figure revealing that an exponent of ~ 1.2 provides 
a fair fit. This is close to the values of 1.10 and 1.35 from 
our relations M bh oc L^°±°- 2 ° and M bh oc L^35±o.30_ 

Thankfully complex models involving the growth of seed 
black holes in high-redshift dark matter halos are not re- 
quired to explain this near-linear relation if it is created 
by the simple addition of similar galaxies. 

4.3. Black hole growth in semi- analytical models 

Given the log-linear relation between Sersic index n 
and spheroid luminosity L (e.g. Graham & Guzman 2003, 
and references therein), the bent Mbh _ fi relation (Gra- 
ham & Driver 2007a) necessitates that the Mbh~£ re- 
lation must, as observed here, also be bent. It would 
appear that black holes within spheroids formed from 
wet-mergers and merger triggered accretion episodes 
(Sanders & Mirabel 1996; Villar-Martm et al. 2012), 
and/or monolithic-like collapse — especially for isolated 
systems — grew via quadratic black hole growth such 
that the Mbh-M sp h C roid relation has a power-law slope of 
2 for the Sersic spheroids. That is, the black hole should 
grow more quickly than the population of stars in the 
host spheroid, and this has in fact recently been observed 
(Seymour et al. 2012; see also Greene et al. 2009). 

Following Kauffmann & Haehnelt (2000), numerous 
theories and models (e.g. Hopkins et al. 2006; Johansson, 
Naab & Burkert 2009; Cen 2012; Hirschmann et al. 2012) 
have been developed to create or match a near-linear 
Mbh-i and Mbh-M sp heroid relation. Unfortunately be- 
cause of this, some of their claimed 'success' implies that 
they are actually in need of refinement. However some 
models may contain important elements or clues that 
have been somewhat over-looked by the community. For 
example, the models of Dubois et al. (2012) reveal that 
while the Mbh~c relation is linear and has a steeper slope 
than 4, their simulated M^-L relation is bent at around 
Mbh = 1O 8 M , as is the simulated Mbh-Ms p heroid rela- 
tion of Cirasuolo et al. (2005) 13 . The hydrodynamical 
simulation models of Khandai et al. (2012, see their fig- 
ure 7) clearly reveal a steep Mbh _ M sp h e roid relation such 

13 Cirasuolo et al. (2005) additionally predict a steepening of the 
M^-cr relation at lower black hole masses. 
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that black hole growth outpaces that of the host spheroid, 
in agreement with Graham (2012a) and this study. 

The primary, so-called 'quasar' or 'cold', mode of black 
hole growth during gas-rich mergers / processes cur- 
rently assumes, in semi-analytical models (Kauffmann & 
Haehnelt 2000, their eq 2; Croton et al. 2006, their eq. 8; 
Guo et al. 2011, their eq. 36) 14 , that black hole growth 
occurs via accretion which is linearly proportional to the 
inflowing mass of cold gas (which also produces the host 
spheroid), modulated by an efficiency which is lower for 
both unequal mass mergers (Croton et al. 2006) and less 
massive (more gas-rich) systems with lower virial veloc- 
ities. We therefore offer the following new prescription 
for the increase in black hole mass, SM^h, due to gas 
accretion during wet mergers: 



SM hh oc 



maj 



l + (280kms ^/Kiriai 



(9) 



Here, M m j n and M ma j are the total baryonic masses from 
the minor and major galaxies involved, and M co id is 
their combined cold gas mass. The term Kiriai is the 
circular or 'virial' velocity of the newly wed galaxy's 
halo, normalized at 280 km s _1 following Kauffmann & 
Haehnelt (2000). This expression, now with an exponent 
of 2 on the cold gas mass, assumes quadratic growth for 
the black hole relative to the stellar mass of the host 
spheroid. 

Another consequence of this order-of- magnitude lower 
Mbh/Mjphcroid mass ratio in faint spheroids is that AGN 
feedback, and in particular the ratio of this feedback 
energy (oc Mbh) to the binding energy of the spheroid 
(oc MgpheroidC 2 ), will be an order of magnitude less than 
previously assumed by studies which used the old con- 
stant value of Mbh/Msphcroid ~ 0.002 (Marconi & Hunt 
2003; Haring & Rix 2004). 

4.4. Barred galaxies 

As noted earlier, barred galaxies can reside up to ~1 
dex below, or perhaps rather to the right of, non-barred 
galaxies in the M^—a diagram, with the mean verti- 
cal offset 0.3 to 0.4 dex (e.g. Graham 2008a; Hu 2008). 
Galaxies with AGN have also been shown to display this 
same behavior in the M^-a diagram (e.g. Graham & 
Li 2009, see also Wandel 2002; Mathur & Grupe 2005; 
Wu 2009; Mathur et al. 2012, while alternate views can 
be found in Botte et al. 2005; Komossa & Xu 2007; De- 
carli et al. 2008a and Marconi et al. 2008). This is not 
to say that all barred galaxies are offset from the bar- 
less Mbh-P relation, only that a significant fraction are. 
In the analysis by Hu (2008), he considered the offset 
galaxies to be 'pseudobulges', although all of his offset 
galaxies were also barred galaxies. Greene, Ho, & Barth 
(2008, see their figure 7b) subsequently showed that their 
sample of 11 'pseudobulges' was not offset from the Mbh~ 
-^spheroid relation defined by massive classical systems 15 . 
(Oddly, their work is often mis-quoted as evidence that 
pseudobulges in disk galaxies are an offset population.) 
With new and updated data, we have observed that the 

14 For clarity we note that Guo et al. (2011) excluded the square 
on the normalized velocity term in their eq. 36. 

15 We speculate that this lack of an offset may have been due to 
an underestimation of dust obscuration from the disks surrounding 
their bulges. 



mean offset of (in this case 21) barred galaxies in the 
Mbh—a diagram remains. These galaxies also appear to 
follow a relation which is parallel with that defined by 
luminous core-Sersic galaxies. The rather narrow base- 
line, of effectively only ^2 mag, for the 11 'pseudobulges' 
in Kormendy et al.'s (2011) M^-L diagram may partly 
(see also Graham 2011b) explain why they did not ob- 
serve a correlation between black hole mass and bulge 
luminosity. It is not yet clear if barred galaxies are offset 
from non-barred (Sersic) galaxies in the Mbh-M sp h e roid 
diagram. 

Graham (2008a) had previously noted the connection 
between bar instabilities and pseudobulges 16 , and wrote 
that if the offset barred galaxies have discrepantly low 
black hole masses rather than elevated cr-values, then 
they should also appear as systematic outliers in the 
M\fo-L diagram. Now although Graham (2007a) 1 , Gra- 
ham & Driver (2007a, their section 3.2) and Graham 
(2008b, his section 2.2.2) warned that the M^-L re- 
lation must be bent 18 , nobody had yet constructed this 
bent relation and it was common practice to assume that 
the original relation, dominated by luminous spheroids, 
could be extrapolated to low luminosities. However in so 
doing, even the classical, low-luminosity spheroids would 
end up incorrectly labeled because they are offset from 
the extrapolated relation fit to the distribution of lumi- 
nous spheroids in the M^-L diagram (Greene, Ho, & 
Barth 2008; Jiang, Greene, & Ho 2011; Kormendy, Ben- 
der, & Cornell 2011). 

In a similar vein, and while on the discussion of pseu- 
dobulge identification, a number of recent papers on the 
topic of supermassive black holes have claimed the de- 
tection of pseudobulges because (due to their low Sersic 
index < 2) they do not follow the extrapolations of other 
linear approximations to the bright arms of what we now 
know are continuous but curved underlying distributions. 
This important subject of non- (log-linear) scaling rela- 
tions is extensively reviewed in Graham (2012c, and ref- 
erences therein) and is vital if we are to properly un- 
derstand the formation of spheroids. For example, semi- 
analytical and theoretical models which still assume a 
log-linear size-luminosity relation (e.g. Guo et al. 2010; 
Suh 2012) are handicapping themselves as they are not 
matching galaxies in the real Universe. Galaxies which 
do not follow the bright arm of the curved (effective sur- 
face brightness)-(effective radius) diagram should also 
not be labeled as pseudobulges for use in the various 
black hole diagrams (e.g. Gadotti & Kauffmann 2009; 
Greene et al. 2008). 

Finally, we offer the following thoughts on whether or 
not barred galaxies might be offset in the M^-L dia- 
gram. Among lenticular disk galaxies, the average lumi- 
nosity of the barred SOs is apparently ^0.4 mag fainter 
than that of the non-barred SOs, largely because there are 
almost no barred SO galaxies brighter than Mb = —21 
mag (van den Bcrgh 2012). It has additionally been re- 

16 Whether or not the bulges in these galaxies formed solely 
from a redistribution of disk stars, to produce 'pseudobulges', is 
unknown due to the many difficulties in identifying pseudobulges 
(Graham 2011a, 2012c). 

17 See the end of his Appendix A, noting that the term L 5 
should have read L 2 . 

18 At this time Bernardi et al. (2007) had speculated that there 
may instead be curvature in the M^-a relation 
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ported that the average photometric bulge mass of barred 
galaxies is less than that of non-barred disc galaxies of 
similar total galaxy stellar mass (Aguerri et al. 2005; 
Laurikainen et al. 2007; Coelho & Gadotti 2011, but 
see Laurikainen et al. 2010), and that the bulge-to-total 
flux ratio is smaller in barred galaxies than in non-barred 
galaxies of the corresponding Hubble type (Weinzirl et al. 
2009). If the effective sizes (R e ) of bulges in barred galax- 
ies are also smaller than those in non-barred disk galax- 
ies, then one may find that the offset of barred galaxies 
toward higher velocity dispersions in the M^-a diagram 
is counter-balanced when using a 2 R c for bulge mass in 
the Mbh--Ms P horoid diagram. Thus, rather than finding 
barred galaxies tending to reside, on average, 'rightward' 
of the Mb^-L relation defined by the non-barred Sersic 
galaxies, the bulges of barred galaxies may yet be found 
to reside 'leftward' of it, i.e. 'above' it. Trying to dis- 
cern if spheroids, and their black holes, in barred galaxies 
are offset from those of non-barred galaxies in the Mbh~ 
L diagram, which was the question posed by Graham 
(2008a) , will require either a larger data sample than we 
currently have or reliable individual bulge/disc/bar de- 
compositions. 

5. CONCLUSIONS 

Using a refined sample of 72 galaxies, comprised of 
24 core-Sersic galaxies plus 48 Sersic galaxies, we have 
revised the Mbh-£ S pheroid relation. We have replaced 
what was a single power-law relation for all spheroids 
with distinct relations for the core-Sersic and the Sersic 
population. 

Using this exact same galaxy sample we have also up- 
dated the Mbh-c relations. From 21 barred and 51 non- 
barred galaxies, we confirm a mean factor of two off- 
set between their M^-a relations, and a steep slope of 
around 5 to 5.5 for the elliptical and un-barred galaxies 
(see Figure EJa, and Table [3}. 

It remains unknown if barred galaxies are offset, to 
lower black hole masses or equivalently brighter lumi- 
nosities, relative to the non-barred Sersic galaxies, in the 



Mbh-ispheroid diagram. 

The lower mass Sersic galaxies are observed to be off- 
set from Mbh _ i S phcroid relation defined by the luminous 
galaxies. However rather than being an offset population 
of alleged pseudobulges, we have explained why it is in- 
stead the luminous core-Sersic galaxies that deviate from 
the Mbh _ i S phcroid relation defined by the fainter Sersic 
galaxies. 

Luminous core-Sersic galaxies roughly follow the rela- 
tion Mbh oc L , consistent with the concept of their for- 
mation through dry galaxy merging. Built by such sim- 
ple additive merger events, rather than what were once 
gaseous mechanisms, a roughly constant (black hole)- 
to-(host spheroid stellar) mass ratio of 0.35-0.55% is 
observed. The fainter (barred and non-barred) Sersic 
galaxies follow an Mbh _ i S pheroid relation which is some 
2.5 times steeper. Thus, while the black hole mass to 
spheroid stellar mass ratio can be 0.5% in massive ellip- 
tical galaxies, it drops to values of 0.05% in spheroids 
that are not much fainter than Mb ~ — 19 mag. 

We have discussed the implications of this steep Sersic 
A^bh-isphcroid relation for the Sersic M bh -Ms p h,dyn re- 
lation (Graham 2012a), which suggest that black holes 
grew quadratically relative to their host spheroids. 
We have therefore advocated a new quadratic cold-gas 
'quasar' mode feeding equation (Eq. [5]) for use in semi- 
analytical models. 

Finally, using our new M b h — -^spheroid relation for Sersic 
galaxies, we have predicted a number of intermediate 
mass black holes, potentially filling a substantial por- 
tion of the current observational gap between stellar mass 
black holes and one million solar mass black holes. 
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